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I.  INTRODUCTION 

A.         BACKGROUND 

In  1 959.  Nonweiler  hypothesized  that  the  pressure  generated  behind  a  planar 
shock  wave  might  be  sufficient  to  provide  the  vehicle  lift.  His  proposal  of  a  caret-shaped 
vehicle  was  the  first  of  a  class  of  vehicles  designed  for  supersonic  flight.  The  principle  of 
lift  generated  by  the  pressure  rise  through  a  shock  wave  lead  to  the  term  "waverider". 
The  decades  that  followed  Nonweiler's  hypothesis  have  brought  a  series  of 
improvements  and  refinements.  The  latest  group  of  proposed  vehicles  is  a  selection  of 
conical  flow  waveriders,  designed  to  capture  and  ride  upon  conical  shaped  shocks.  A 
number  of  studies  on  the  optimization  of  waverider-configured  vehicles  are  ongoing 
[Refs.  1  and  2]. 

The  Naval  Postgraduate  School  (NPS)  is  pursuing  in  an  ongoing  conceptual 
development  of  hypersonic,  carrier-based  waverider-configured  aircraft.  Since  1992 
both  individual  thesis  students  and  design  team  efforts  have  contributed  to  the  body  of 
knowledge  in  this  emerging  field  of  aerodynamics  [Ref.  3].  Included  within  this  study  is 
the  development  and  investigation  of  a  proposed  carrier-based  interceptor  aircraft,  and 
among  the  proposed  designs  is  the  Price  waverider. 

Using  the  NASA  Ames  Research  Center  (ARC)  Waverider  Code,  David  Price 
developed  a  conceptual  configuration  for  a  hypersonic,  carrier-based,  deck-launched, 


waverider  strike  aircraft,  optimized  for  Mach=6, 


D^ 


,  flight  [Ref.  4].  The 


waverider  code  is  a  subset  of  the  Hypersonic  Aircraft  Vehicle  Optimization  Code 
(HAVOC),  created  by  the  Systems  Analysis  Branch  at  NASA  ARC.  Optimization  of  the 
product  of  L/D  and  the  specific  impulse  produced  an  integrated  engine  and  air  vehicle. 
This  method  is  in  concert  with  one  of  the  fundamental  design  aspects  of  hypersonic 
vehicles.  The  first  and  most  important  aspect  is  to  create  and  capture  a  conical  shock.  A 
second  important  aspect,  however,  capitalizes  on  shock  compression,  forward  of  the 
engine  intake.  This  shock  pre-compresses  the  engine  inlet  air.  Aft  body  shaping  is  used 


to  expand  the  exhaust  in  an  external  nozzle.  Thus  engine  air  compression  and  expansion 
takes  place  external  to  the  engine  modules  themselves  allowing  for  a  reduction  in  weight, 
complexity  and  cooling  requirements.  In  this  sense  the  engine  and  aircraft  configuration 
must  be  an  integrated  effort,  both  being  interdependent  in  all  aspects  of  flow  analysis  - 
hence  the  need  to  optimize  this  type  of  vehicle  using  both  aerodynamic  parameters  (L/D) 
and  engine  parameters  (Isp). 

Follow-on  work  to  Price's  initial  development  are  ongoing.  Based  on  this 
configuration,  15-inch  and  8-inch  root  chord  aluminum  models  were  produced  by  NASA 
Ames  for  NPS  research  to  facilitate  the  further  characterization  of  the  performance  of  the 
Price  waverider.  Using  the  smaller  aluminum  model,  Lowell  M.  Johnson  performed  a 
series  of  flow  visualization  tests  of  the  Price  waverider  in  the  NPS  water  tunnel  [Ref.  5]. 
Subsonic  wind  tunnel  testing  of  the  Price  waverider  has  been  performed  by  Mark  E. 
Cedrun  [Ref.  6]  and  Michael  Huff  [Ref.  7],  using  the  larger  of  the  aluminum  models. 
Tests  to  date  indicate  that  flight  is  feasible  for  waverider-conflgured  aircraft.  The  Price 
waverider  database,  while  including  water  and  wind  tunnel  data,  lacks  a  Computational 
Fluid  Dynamic  (CFD)  component.  Additionally,  no  supersonic  or  hypersonic  data  were 
available  for  the  Price  waverider  prior  to  this  CFD  analysis. 

The  primary  thrust  of  this  thesis  was  the  development  of  the  appropriate  surface 
grid  for  the  Price  waverider  and  the  calculation  of  aerodynamic  characteristics  in 
subsonic  and  supersonic  flight.  The  subsonic  CFD  flow  solutions  were  aimed  at 
complimenting  existing  wind  tunnel  and  water  tunnel  data.  CFD  results  are  compared 
with  the  force  and  moment  data  obtained  in  recent  wind  tunnel  tests  for  a  range  of 
Reynolds  numbers  and  angles  of  attack.  The  range  of  subsonic,  transonic  and  supersonic 
solutions  represents  a  continuation  of  the  ongoing  research  effort  supporting  the  complete 
performance  analysis  of  the  Price  waverider. 

B.         OBJECTIVE 

The  objectives  of  this  thesis  are  to  determine  the  theoretical  force  and  moment 
data  over  a  broad  Mach  number  spectrum,  and  to  compare  theoretical  and  experimental 
results  in  the  subsonic  flight  regime.  The  results,  across  the  whole  flight  spectrum  of 


Mach  numbers,  are  to  serve  as  a  basis  for  continued  examination  and  feasibility  analysis 
of  the  Price  waverider. 

C.         PRICE  WAVERIDER 

The  Price  waverider  was  optimized  for  Mach  6.0  flight  using  inviscid,  conical 
shock,  optimization  methods  contained  in  the  Ames  Research  Center  Waverider  Code. 
The  parabolic  leading  edge  droops  from  the  nose  to  the  trailing  edge.  While  the  upper 
surface  remains  entirely  in  line  with  the  freestream,  the  lower  surface  slopes  to  create  a 
wedged-body  over  which  the  shocks  are  formed.  Along  the  centerline,  a  pair  of  ramps 
precedes  the  engine  inlet.  Off  the  centerline  the  lower  surface  creates  a  5.2  degree  semi- 
vertex  angle.  Along  the  aft-body  the  lower  surface  is  faired  upward  to  meet  the  trailing 
edge.  Waveriders  derived  from  strictly  conical  flow  will  have  a  flat  base  at  the  trailing 
edge,  creating  enormous  base  drag,  thus  the  lower  surface  must  be  faired  to  narrow  the 
trailing  edge  to  decrease  the  base  drag  or  pressure  drag  of  the  vehicle  [Ref.  8].  The  fairing 
of  the  aft  body  under-surface  up  to  the  trailing  edge  creates  a  negative  camber  in  the 
body.  This  negative  camber  is  more  pronounced  on  the  centerline,  where  the  engine 
cowling  ramps  back  up  to  the  trailing  edge. 
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n.         DESCRIPTION  OF  CFD 

Computational  Fluid  Dynamic  methods  were  developed  as  a  means  of  solving  the 
complex  system  of  partial  differential  equations  that  describe  a  given  flowfield.  The  use 
of  computers  to  solve  discretized  equations  enables  the  solution  to  be  found  for  those 
partial  differential  equations  for  which  no  exact  solution  exists.  Since  1975,  numerous 
advancements  have  been  made  to  refine  the  application  of  several  finite  element  methods 
for  computing  flowfields;  improving  algorithms  to  improve  computational  efficiency  and 
further  defining  and  extending  the  regimes  in  which  CFD  methods  can  be  applied  with 
confidence  [Ref.  3]. 

The  application  of  CFD  techniques  to  a  known  geometry  is  a  multi-step  process. 
Several  software  applications  are  used  for  grid  generation  and  refinement,  and  flow 
solution  and  analysis.  NPS  has  a  suite  of  computer  software  that  enables  grid  generation 
and  manipulation,  flow  solving  and  graphic  presentation  of  the  solution. 

A.  GRIDGEN  OVERVIEW 

GRIDGEN  Version  9  is  a  Computer  Science  Corporation  developed  software 
tool,  which  was  originally  sponsored  for  development  by  the  U.  S.  Air  Force.  It  is  an 
interactive  code  used  to  create  2-,  or  3-dimensional  grids.  The  domains  can  be 
represented  graphically  within  GRIDGEN  or  exported  as  data  files  for  additional 
operations,  such  as  providing  the  computational  space  for  flow  solutions.  GRIDGEN 
requires  the  import  of  surface  points  from  a  previous  design  effort  or  the  creation  of 
surfaces  manually  within  the  bounds  of  the  program  [Ref.  9]. 

B.  GRIDGEN3D  OVERVIEW 

GRIDGEN3D  is  a  batch  code  used  in  conjunction  with  GRIDGEN  Version  9. 
This  code  is  used  to  refine  the  grid  volume  by  adjusting  grid  points  internal  to  the  block 
of  points.  In  using  of  GRIDGEN3D,  caution  should  be  applied  to  ensure  that  the  body 


points  are  fixed  to  prevent  the  blending  of  body  features  by  designating  body  surface  as 
fixed  boundaries. 

C.  GRIDED 

GRIDED  is  an  interactive,  menu-driven  file  FORTRAN-based  software 
application  used  to  manipulate  the  grid.  GRIDED  allows  the  addition  and  removal  of 
grid  points  or  planes,  the  selection  of  indices  and  the  performance  of  coordinate 
transforms. 

D.  OVERFLOW1.8B  OVERVIEW 

OVERFLOW  1.8b  is  a  flow  solution  code  developed  by  NASA  Ames  Research 
Center  and  is  a  combination  of  the  older  F3D  and  Chimera  codes.  It  uses  an  implicit, 
finite-differencing  scheme  to  solve  the  Reynolds-Averaged  Navier-Stokes  equations  in 
strong  conservative  form  [Ref.  10].  The  input  file  from  the  user  enables  the  selection  of: 

•  Flight  regime  or  flow  conditions,  such  as  Mach  number,  angle  of  attack,  yaw 
angle,  ratio  of  specific  heats,  viscosity  and  freestream  temperature. 

•  Turbulence  model  selection  (Options  include  Baldwin-Lomax,  Baldwin- 
Barth,  k-G3  and  k-e  models) 

•  Computational  controls  such  as  time  step,  differencing  scheme,  artificial 
viscosity  selection  and  smoothing  parameters. 

•  Boundary  conditions  (options  include  viscous  or  inviscid  walls,  inflow  or 
outflow,  or  flow-through  cuts). 

OVERFLOW  1.8b  solutions  are  run  on  the  Cray  Y-MP  J94  computer,  owned  and 
operated  by  the  Naval  Postgraduate  School. 

E.  FAST 

Flow  Analysis  Software  Toolkit  (FAST)  was  used  to  compute  and  graphically 
represent  the  grid  and  its  corresponding  flow  solution.  The  images  of  the  waverider  body 
and  the  Mach  number  and  pressure  profiles  contained  within  Chapter  IV  are  products  of 
FAST.  [Ref.  1 1  ] 


ni.       GRID  GENERATION  AND  FILE  INPUT 

Development  of  an  appropriate  grid  and  the  correspondingly  appropriate 
boundary  conditions  and  flowfield  constraints  was  essential  to  the  production  of  realistic 
solutions  for  the  waverider  flowfield.  Additionally,  a  number  of  options  exist  within 
OVERFLOW1.8b,  the  flow  solver,  for  refining  the  techniques  used  to  optimize  the 
computational  efficiency  or  to  enhance  certain  flow  characteristics. 

A.         GRID  GENERATION 

Grid  generation  was  accomplished  using  the  previously  described  code, 
GRIDGEN.  Using  grid  points,  provided  by  NASA  Ames  (originally  generated  by  David 
Price),  the  surface  geometry  of  the  Price  Waverider  was  defined  by  importing  the  series 
of  upper-surface  and  lower-surface  fuselage  station  grid  points  into  GRIDGEN.  The 
surface  domains  were  then  built  up  by  manually  inserting  the  coordinates  in  GRIDGEN 
to  create  sequential  fuselage  station  connectors.  Once  both  the  upper  and  lower  surface 
of  the  waverider  was  completely  defined  by  the  series  of  linked  fuselage  station  domains, 
these  domains  were  consolidated  into  a  single-domain  covering  the  upper  and  lower 
surfaces.  Following  the  initial  creation  of  the  waverider 's  surface  grid,  a  number  of 
anomalies  were  observed  in  GRIDGEN  in  the  basic  geometry  of  the  grid,  so  minor 
adjustments  to  the  database  provided  by  NASA  were  made. 

Those  changes  included  a  coordinate  transformation,  elimination  of  co-incident  or 
nearly  co-incident  fuselage  stations  and  minor  body  modifications  to  match  model 
geometry.  The  fuselage  station  coordinates  that  were  provided  by  NASA  Ames  created  a 
body  whose  axis  lay  in  the  z-direction.  Since  OVERFLOW1 .8b  specified  Mo  in  the  x- 
direction,  a  coordinate  transformation  was  performed  to  a  right-handed  coordinate  system 
with  the  axis  of  the  body  lying  on  the  x-axis.  After  an  examination  to  ensure  there  was 
no  loss  of  geometric  features  on  the  waverider  surface,  fuselage  stations  at  x  =  0.03014, 
0.03024,  0.03034,  0.33192,  0.63181  and  0.65947,  were  eliminated  due  to  their  positions, 
nearly  coincident  to  adjacent  stations.  Additionally,  the  grid  point  defining  the  outermost 
point  of  the  fuselage  station  at  x  =  0.93164  appeared  to  be  missing,  creating  a  sawtooth  in 


the  leading  edge.  A  new  end  point  was  created  for  this  fuselage  station  to  bring  the 
points  at  that  station  out  to  create  a  smooth  leading  edge,  rather  than  the  sawtooth 
observed  in  GRIDGEN.  The  new  point  was  created  by  interpolating  between  the  outer 
points  of  adjacent  fuselage  stations.  Lastly,  the  engine  exhaust  in  the  wind  tunnel  model 
was  a  ramped  exhaust,  whereas  the  points  provided  by  NASA  Ames  modeled  the  exhaust 
as  a  flat  plate.  To  match  the  geometry  of  the  wind  tunnel  model  and  to  match  a  more 
realistic  full-scale  configuration,  the  fuselage  station  at  x  =  0.93816  was  eliminated  in  its 
entirety.  The  resulting  engine  exhaust  ramp  geometry  then  mimicked  the  model.  The 
surface  grid  points,  with  the  required  adjustments,  are  provided  in  Appendix  A. 

The  spanwise  grid  dimension  (grid  dimensions  refer  to  the  number  of  points  in  the 
grid  in  a  given  direction)  of  the  waverider  surface  was  set  at  20  grid  points  each  for  the 
upper  and  lower  surfaces.  After  some  experimentation,  20  grid  points  in  the  span  was 
seen  as  the  minimum  number  of  points  required  to  model  the  surface  geometry  clearly. 
Additional  points  would  increase  the  size  of  the  grid  and  the  grid  density,  slowing 
computations.  Clustering  points  at  the  outer  edge  of  the  body's  upper  surface  was  done 
to  provide  sufficient  grid  point  density  near  the  sharp  leading  edge.  The  body's  lower 
surface  required  grid  point  clustering  at  both  the  leading  edge  and  at  the  engine  cowling 
in  order  to  ensure  adequate  modeling  of  the  body  geometry.  The  chord  wise  dimension  of 
the  grid  was  32  points  from  nose  to  tail,  after  the  removal  of  the  fuselage  stations  noted 
above.  Additional  gridlines  were  added  to  provide  higher  grid  density  at  the  trailing  edge 
of  the  engine  cowling  and  at  the  trailing  edge  of  the  vehicle. 

After  numerous  attempts  to  eliminate  negative  cell  volumes  at  the  trailing  edge,  a 
cusp  was  added  at  the  trailing  edge  to  ease  grid  geometry.  This  addition  extended  the 
blunt  trailing  edge  a  distance  of  1 .5  percent  of  the  root  chord  length.  This  cusp  is  the 
only  modification  made  that  changed  the  computational  grid  shape  from  the  model  shape. 
Creation  of  a  sharply-cusped  trailing  edge  was  believed  to  be  more  representative  of  final 
vehicle  shape,  as  this  addition  helped  to  alleviate  the  high  base-pressure  drag  forces 
caused  by  the  bluntness  of  the  Price  waverider  trailing  edge.  The  final  dimension  of  the 
grid  was  45  points  from  nose  to  tail.  Overall  length  of  the  vehicle,  as  prescribed  by  the 


grid  points  provided  by  NASA  Ames,  was  equal  to  1.0,  prior  the  addition  of  the  cusped 
trailing  edge.  Dimensions,  when  quoted  as  a  fraction  of  the  body  length,  are  a  fraction  of 
the  original  body  length  or  root  chord,  vice  a  fraction  of  the  computational  model  length 
or  root  chord. 

The  resulting  surface  grid  is  shown  in  Figures  1  and  2.  The  inset  of  Figure  1 
shows  the  addition  of  the  trailing  edge  cusp,  increasing  the  total  model  length  by  1.5 
percent. 


dx  =  0.015 


Figure  1 .  Waverider  Grid  Three-View  and  Oblique 


Figure  2.  Waverider  Surface  Grid  Cusp 

In  addition  to  the  complex  geometry  of  the  lower  surface  and  the  sharp  curvature 
of  the  trailing  edge  cusp,  the  droop  of  the  leading  edge  of  the  wing  contributed  to  the 
difficulty  in  creating  a  single  three-dimensional  block  that  was  suitable  for  use  by 
OVERFLO  Wl  .8b.  The  vehicle  configuration  required  a  multi-step  approach  to  creating 
a  suitable  3 -dimensional  computational  block  surrounding  the  waverider  half-body.  First, 
it  was  necessary  to  split  the  vehicle's  upper  and  lower  surfaces  and  map  out  the 
surrounding  airspace  in  two  separate  blocks  to  prevent  gridlines  from  crossing  from  the 
lower  surface  to  the  upper  surface.  Then  the  two  blocks  had  to  be  merged  to  create  the 
entire  computational  space.  This  was  done  by  creating  the  surrounding  surfaces 
including  domains  with  the  following  number  of  grid  points: 


Sting: 

20x50 

Wake  Cut: 

20x50 

Cut  Plane: 

45x50 

Outer  Upper  Shell: 

20x45 

Outer  Lower  Shell: 

20x45 

Upper  Back  Plane: 

45x50 

Lower  Back  Plane: 

45x50 

Of  those  domains,  the  sting,  the  wake,  and  the  plane  that  was  created  to  divide  the 
upper  and  lower  surface  grids  were  included,  identically,  in  both  the  upper  and  lower 
blocks  to  be  joined  as  identical  elements.  It  should  be  noted  that,  in  this  application,  the 
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sting  is  the  representation  of  the  line  forward  of  the  waverider  nose,  created  as  an  entity 
of  the  computational  space.  It  does  not  refer  to  the  sting  used  in  wind  tunnel  tests  to 
mount  the  vehicle  in  the  tunnel  test  section.  From  these  domains,  along  with  the  upper 
and  lower  surfaces  of  the  body,  the  six-domain  upper  and  lower  3D  volumetric  grids, 
known  as  blocks,  were  created.  Inspection  of  each  of  the  blocks  using  GRIDGEN3D, 
was  used  to  determine  the  presence  of  negative  cell  volumes.  Negative  cell  volumes 
resulted  when  the  grid  lines  within  the  grid  volume  had  tangled  and  the  indices  of  the 
grid  lines  in  certain  areas  were  no  longer  sequential.  The  implication  of  this  is  that  Ax  or 
Ay  used  by  the  finite  differencing  scheme  in  OVERFLOW  1.8b  will  be  a  negative 
number,  which  is  unacceptable  to  the  code. 

The  GRIDGEN3D  elliptic  solver,  with  a  LaPlacian  foreground  control  function, 
was  used  to  resolve  the  negative  cell  volumes  in  both  the  upper  and  lower  blocks.  The 
LaPlacian  function  tended  to  drive  the  cells  within  the  3D  block  to  equal  volume  cells 
throughout  the  total  computational  volume,  thus  those  cells  with  negative  volumes  were 
driven  to  match  the  average  cell  volume,  which  was  positive.  While  equal  volume  was 
unlikely  to  be  achieved  in  a  3D  grid  containing  a  sharp  body,  applying  this  control 
eliminated  the  crossing  of  gridlines.  It  was  possible  to  allow  the  elliptic  solver  to  run  for 
too  many  iterations.  If  this  was  done  the  number  of  negative  cell  volumes  tended  to 
increase.  Periodically  the  elliptic  solver  was  stopped,  and  the  grid  was  examined  to 
ensure  the  number  of  cells  with  negative  volume  were  continuing  to  decrease. 
Regardless  of  the  care  taken  not  to  over  manipulate  the  grid,  the  elliptic  solver  tended  to 
increase  the  distance  between  the  body  and  its  adjacent  grid  points  at  all  points  except  the 
plane  of  symmetry,  where  grid  points  were  held  as  fixed  entities. 

FAST  images  of  the  upper  and  lower  3D  blocks,  following  the  elimination  of  the 
negative  cell  volumes,  were  examined  to  identify  coordinate  discrepancies.  It  was 
necessary  to  reverse  the  direction  of  the  J-planes  (the  surfaces  radiating  from  the  body  to 
the  outer  shell  in  lengthwise  planes)  and  K-planes  (surface  planes  stepping  out  from  the 
body  surface  to  the  outer  shell)  of  the  upper  block.  This  was  accomplished  using  the 
GRIDED  option  to  reverse  plane  indices.  Following  coordination  of  the  plane  indices  of 
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the  upper  and  lower  blocks,  the  two  blocks  were  joined,  again  using  GRIDED. 
Additionally,  one  J-plane  was  added  to  the  upper  and  lower  back  planes  of  the  block  to 
provide  the  planes  of  symmetry  that  facilitated  the  use  of  a  central-differencing  scheme  in 
the  computational  phase  of  this  analysis.  This  single  block  was  then  saved  and  examined, 
once  again,  for  negative  cell  volumes  using  GRIDGEN3D. 

The  figure  below  provides  a  perspective  view  of  the  3-D  computational  block.  At 
the  center  of  this  block  is  the  surface  of  the  waverider,  seen  as  the  black  body.    Final 
dimensions  of  the  computational  block  were  41x45x50  grid  points,  for  a  total  of  94,250 
points. 


Figure  3.  Three-Dimensional  Computational  Block 
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Coordinate  directions  used  by  OVERFLOW1.8b  were  specified  to  have  the  J- 
direction  increasing  along  the  length  of  the  body  from  nose  to  tail.  The  K-direction  was 
specified  as  positive  around  the  body,  from  the  upper  centerline,  out  to  the  wing  tip,  and 
back  into  the  centerline  of  the  lower  surface.  The  L-direction  is  positive  in  the  direction 
radiating  out  from  the  body  to  the  shell  of  the  3D  grid. 

B.         SELECTED  FLIGHT  CONDITIONS 

To  support  the  study  of  Reynolds  number  effects  for  both  viscous  and  inviscid 
flow  conditions,  a  range  (0.3  <  Moo  <  6.0)  of  computations  were  conducted.  A  summary 
of  the  flight  conditions  for  which  analysis  was  performed  is  provided  in  Table  1. 


Mach  Number  Range 

Angle  of  Attack 
Sweep 

Engine  Flow- 
Through 

0°AOA 

Mo  =  0.3 

0°AOA 

Inviscid 

Viscous 

Viscous 

Inviscid 

0.3 

0.3 

(degrees  AOA) 

Mo  =  6.0 

0.4 

0.4 

0.5 

0.5 

-10 

0.6 

0.6 

-8 

0.7 

0.7 

-4 

0.8 

0.8 

-2 

0.9 

0.9 

0 

1.2 

1.2 

6 

2.0 

2.0 

2 

4.0 

4.0 

8 

6.0 

6.0 

10 

Table  1.  Summary  of  Flight  Conditions  Analyzed 

Both  viscous  and  inviscid  flow  solutions  for  the  Price  waverider  were  run  for  Moo 
=  0.3,  0.4,  0.5,  0.6,  0.7,  0.8,  0.9,  1.2,  2.0,  4.0,  and  6.0.  Alpha  sweeps  were  conducted  at 
Moo  =  0.3.    Additionally,  the  Moo  =  6.0  input  file  was  modified  to  include  boundary 
conditions  simulating  engine  intake  and  exhaust  features  to  examine  the  effect  of  engine 
flow-through  at  design  conditions.  OVERFLOWl  .8b  is  optimized  for  Moo  -  5.0  flow 
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solutions.  While  the  code  performed  well  over  a  broad  range  of  Mach  numbers,  the  time 
required  to  compute  low  speed  solutions  and  the  decreasing  stability  of  the  code  at  Mach 
numbers  greater  than  Moo  =  5.0  began  to  limit  its  application. 

Figure  4  shows  the  number  of  iterations  required  to  achieve  convergence  in  the 
flow  solutions  at  each  Mach  number.  At  the  lower  Mach  number,  and  indeed  even  at  M„ 
=  0.3,  the  computational  time  required  to  fully  converge  a  solution  is  over  130,000 
iterations.  This  equated  to  over  50  hours  of  computational  time,  using  all  four  processors 
of  the  Cray  supercomputer.  The  Mach  0.5  solution  required  nearly  160,000  iterations  due 
to  the  smaller  time  step  that  was  required  to  keep  the  solution  stable.  The  reason  for  the 
instability  seen  at  Moo  =  0.5  was  not  clear.  At  M„  =  6.0,  the  time  step  required  to  prevent 
the  solution  from  becoming  unstable  was  very  small  and  the  computational  solution  was 
particularly  intolerant  of  even  moderately  aggressive  attempts  to  speed  up  the 
convergence. 
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Figure  4.  Number  of  Iterations  vs.  Mach  Number 
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C.         FLOW  SOLUTION  BATCH-FILE  INPUT 

OVERFLOW  1.8b  is  a  batch-run  FORTRAN-based  flow  solver,  which  required 
the  creation  of  input  files  to  start  each  flow  solution.  Input  files  for  each  flight  regime 
were  created,  using  Reference  8,  to  reflect  the  flow  conditions,  desired  turbulence  model, 
computational  controls  and  boundary  conditions.  Each  of  the  batch  input  files  was  set  up 
to  apply  a  central-differencing  scheme  on  a  single  grid  computational  space.  Mach 
number  and  Reynolds  number  inputs  reflected  the  desired  freestream  Mach  number  and 
the  Reynolds  number  computed  for  the  15-inch  model  at  sea-level  conditions.  Boundary 
conditions  at  the  upper  and  lower  back  planes  of  the  3D  grid  were  defined  as  planes  of 
symmetry  in  the  Y-direction.  The  wake  cut  was  intended  to  be  invisible  to  the  flow  and 
was  defined  as  a  fold-over  cut  flow-through  plane.  The  outer  shell  was  defined  as  a 
supersonic  and/or  subsonic  inflow/outflow  surface. 

Throughout  the  course  of  computations  the  input  files  were  updated  to  reflect 
more  aggressive  time  steps,  allowable  after  a  solution  began  to  converge.  The  time  steps 
were  adjusted  without  regard  to  real-time.  Sample  input  batch-files  are  provided  in 
Appendix  B. 

A  summary  of  the  general  approaches  used  is  given  below. 

1.         Subsonic  versus  Supersonic  Flowfields 

In  general,  the  allowable  time  steps  for  the  finite  difference  calculations  were  a 
function  of  the  freestream  Mach  number  and  maturity  (i.  e.,  how  far  the  solution  had 
converged)  of  the  solution.  Once  the  solution  began  to  converge  time  steps  could  be 
increased,  usually  by  factors  of  two  until  convergence.  Allowable  time  steps  were 
greatly  reduced  at  M,x,  =  6.0  because  of  the  velocity  of  the  flow.  The  time  step  here  had 
to  be  reduced  by  an  order  of  magnitude,  to  dt  =  0.001,  to  meet  stability  requirements. 
Tannehill,  Anderson  and  Pletcher  provide  an  excellent  discussion  of  stability 
requirements  and  Courant-Friedrich-Levy  (CFL)  numbers  [Ref.  12] . 

Transonic  and  supersonic  flow  solutions  required  the  use  of  second-order 
smoothing  in  all  three  directions  to  accommodate  discontinuities  in  the  flowfield  across 
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shocks.   Shock  discontinuities  became  present  in  the  flowfield  at  freestream  Mach 
numbers  of  Moo  =  0.7  and  higher.  Prior  to  that  Mach  number,  no  special  considerations  in 
the  input  files  were  required. 

2.  Inviscid  versus  Viscous  Flowfields 

The  input  files  for  Euler  (inviscid)  solutions  required  fewer  arguments  as  the  flow 
description  defined  no  viscosity  effects  in  any  direction  nor  was  any  turbulence  model 
specified.  Conversely  the  input  batch-files  for  the  viscous  cases  specified  viscosity 
effects  in  the  J-  (around  the  waverider  body)  and  L-directions  (from  the  waverider  body 
to  the  outer  shell).  Viscosity  in  the  K  direction  (along  the  length  of  the  waverider  body) 
was  neglected  in  order  to  reduce  the  computational  time.  Neglecting  these  streamwise 
viscous  terms  resulted  in  the  solution  of  a  parabolized  form  of  the  Navier-Stokes 
equations,  which  was  computationally  more  efficient.  The  Baldwin-Lomax  turbulence 
model  was  selected  for  use  in  all  viscous  cases. 

3.  Engine  Flow-Through  Modeling 

Input  batch-files  for  the  engine  flow-through  for  required  additional  boundary 
conditions  that  described  the  intake  area  as  a  region  of  extrapolated  outflow  (i.  e.,  the 
flow  approaching  the  intake  was  allowed  to  exit  the  computational  domain  through  that 
region,  shown  in  Figure  5).  This  figure  shows  the  flow  approaching  the  engine  inlet  from 
the  right  and  penetrating  the  area  of  the  inlet,  seen  as  arrows  inside  the  engine  duct. 
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Figure  5.  Engine  Flow-Through  Model  Intake 

The  engine  exhaust  are  was  defined  as  nozzle  inflow  (i.  e.,  the  mass  flow, 
extrapolated  from  po  and  To,  was  injected  into  the  computational  domain  at  the  exhaust 
plane).  Engine  exhaust  flow  was  defined  as  freestream  flow,  shown  in  Figure  6.  Here 
we  see  both  the  results  of  the  external  flow  solution  as  green  arrows  and  the 
superimposed  engine  flow-through  exhaust,  seen  as  the  arrows  in  the  freestream 
direction. 


Figure  6.  Engine  Flow-through  Model  Exhaust 
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IV.       RESULTS  AND  DISCUSSION 

A.         SUBSONIC  RESULTS 

Integral  to  the  validation  of  computational  solutions  is  the  attainment  of  "steady- 
state"  or  convergence.  As  the  boundary  conditions  propagate  through  the  computational 
space,  values  at  a  given  point  are  compared  to  the  values  of  the  previous  iteration  until 
the  residual  difference  becomes  negligible.  In  this  study,  a  drop  in  residual  values  of 
three  orders  of  magnitude  was  considered  a  converged  solution.  A  sample  convergence 
history  is  provided  in  Figure  7. 

To  validate  each  of  the  flow  solutions  reported  in  this  study,  each  flight  condition 
solution  was  examined  for  convergence  of  the  solution's  residuals,  convergence  of  the  lift 
and  drag  coefficients  to  a  steady  state,  and  representative  surface  pressure  distribution 
and  Mach  number  distribution  within  the  flowfield  surrounding  the  vehicle. 

1.  Mach  0.3  In  vise  id  Results 

The  significant  computational  time  required  to  conduct  low  subsonic 
computational  analysis  limited  low-subsonic  solutions  to  a  minimum  Mach  number  of 
Moo  =  0.3.  Wind  tunnel  tests  provided  data  up  to  a  maximum  Mach  number  of  Moo  = 
0.15.  While  the  computational  analysis  provided  here  is  not  entirely  coincidental  with 
existing  wind  tunnel  data,  the  solution  obtained  demonstrated  fair  representation  of  the 
low  speed  flow  characteristics  of  the  Price  waverider. 

Examining  the  Moo  =  0.3  inviscid  solution,  shown  in  Figure  7,  shows  that  the 
convergence  history  achieved  just  over  three  orders  of  magnitude  convergence  in  126,200 
iterations.  The  short  period  oscillations  were  the  result  of  periodic  vortex  formation  and 
subsequent  shedding  at  the  trailing  edge  of  the  vehicle.  Also  present  in  the  convergence 
history  was  a  long  period  oscillation  occurring  approximately  once  every  thousand 
iterations.  This  was  speculated  to  be  the  result  of  acoustic  waves  reverberating  within  the 
computational  space.  While  no  direct  evidence  within  the  solution  of  the  flowfield 
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supported  that  conclusion,  their  presence  in  the  subsonic  solutions  only  supported  the 
conclusions  that  an  acoustic  phenomena  of  some  type  had  occurred. 
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Figure  7.  Mach  0.3  (Inviscid)  Residual  Convergence  History 

Figure  8  shows  the  drag  coefficient  history.  Note  that  the  drag  coefficient  also 
converged  to  a  steady  state  as  the  solution  converged.  The  final  achieved  drag  coefficient 
was  computed  was  Co  =  0.0194. 


20 


0.18 


0.16  - 


0.14  - 


0.12   - 


0.1 


0.08 


0.05- 


0  04  - 


0.02  - 


1^  =  0.3 
I n viscid  Solution 


1 


CD  =  0.0194 


6  8 

Hecalions 


10 


12 


zlO 


Figure  8.  Drag  Coefficient  vs.  Number  of  Iterations  for  M^  =  0.3,  Inviscid 

Figure  9  shows  the  lift  coefficient  history.  This  force  coefficient  converged  to  a 
stable  value  coincident  with  the  drop  in  residuals.  At  M,  =  0.3,  the  lift  coefficient 
converged  to  Cl  =  -0.03 10.  The  presence  of  a  negative  lift  coefficient  at  subsonic 
speeds  was  the  result  of  the  body  geometry  that  produced  lift  at  supersonic  speeds.  While 
the  wedge  angle  of  the  vehicle  forebody  created  a  mean  chord  line,  that  was  a  straight  line 
at  a  positive  angle  of  attack,  the  nozzle  expansion  aft  of  the  engine  cowl  and  the  shape  of 
the  lower  wing  surface  created  a  negative  camber  in  the  vehicle,  giving  rise  to  negative 
lift  at  zero  degrees  angle  of  attack. 


21 


02 


0.15 


0.05  - 


-O.05  " 


-0.1  - 


-0.15  - 


-02  - 


-0.25 


Moo  =  0.3 
Inviscid  Solution 


CL  =  -0.031 


6  8 

leralions 


10 


12 


xiO 


Figure  9.  Lift  Coefficient  vs.  Number  of  Iterations  for  Mo  =  0.3,  Inviscid 

Using  FAST  images,  the  flowfield  solution  was  graphically  examined.  Figure  10 
shows  the  distribution  of  the  normalized  pressure  (pressure  is  normalized  to  the 
freestream  static  pressure)  on  the  body  surface  and  the  Mach  number  contours  normal  to 
the  body  in  the  surrounding  flowfield.  The  normalized  pressure  was  equal  to  the 
freestream  pressure  over  the  entire  upper  surface.  In  the  inviscid  case,  the  upper  surface 
flow  was  unperturbed,  as  the  waverider  upper  surface  was  aligned  with  the  freestream. 
The  lower  surface  pressure  distribution  illustrated  a  slightly  higher  pressure  area  on  the 
lower  surface  of  the  vehicle  forebody  than  on  the  lower  aftbody,  indicating  a  potential  for 
significant  pitch-up  in  the  vehicle  at  low  speeds.  This  was  borne  out  by  the  computed 
pitching  moment  coefficient  of  Cm  =  -0.033.  Neither  region  displayed  a  pressure 
distribution  that  was  capable  of  providing  lift,  as  was  confirmed  by  the  negative  lift 
coefficient. 
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Figure  10.  Normalized  Surface  Pressure  and  Mach  Number  Contours,  Mo  =  0.3,  Inviscid 

The  Mach  number  contours  confirmed  the  nearly  unperturbed  upper  surface, 
which  contrasted  with  the  large  Mach  number  gradients  at  the  trailing  edge  and  the  less 
extreme  but  critical,  Mach  number  contours  along  the  lower  surface.  The  Mach  number 
contours  at  the  nose  of  the  vehicle  indicated  that  the  stagnation  point  was  located  on  the 
lower  surface  of  the  body  slightly  aft  of  the  nose  of  the  vehicle.  Contours  along  the  lower 
surface  outlined  the  flow  compression  at  the  beginning  of  the  first  and  second  engine 
inlet  ramps  and  subsequent  expansion  as  the  flow  turned  back  toward  the  freestream  at 
the  engine  cowl  and  upward  again  at  the  engine  nozzle.  Further  away  from  the  body,  grid 
spacing  was  less  dense,  which  resulted  in  the  Mach  contour  lines  not  being  represented  as 
smooth  curves.  This  lack  of  smooth  contour  lines,  indicated  a  grid  dependency  within  the 
solution. 
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2.         Mach  0.3  Viscous  Results 

Viscous  terms  were  added  to  the  partially  converged  inviscid  solution  at  66,200 
iterations  in  order  to  capitalize  on  the  convergence  already  achieved  during  the  inviscid 
computation.  OVERFLOW  1.8b  allowed  the  viscous  modification  to  previously  inviscid 
solutions.  This  addition  had  several  notable  impacts.  Shown  in  Figure  11,  residuals  rose 
sharply  as  the  flow  was  "virtually"  perturbed  by  the  addition  of  viscosity,  spiking  at  well 
over  the  initial  maximum,  but  recovering  within  1000  iterations.    As  is  seen  in  Figure  11, 
short  period  oscillations  remained  present,  as  the  more  physically  representative  wake  in 
the  viscous  model  continued  to  generate  and  shed  vortices.  The  sonic  wave  reverberation 
persisted  as  well.  It  should  be  noted  that  the  addition  of  viscous  terms  enabled  a  slightly 
more  aggressive  time  step  to  be  used  in  the  remaining  computations. 

It  was  hoped  that  the  addition  of  these  viscous  terms  would  serve  to  stabilize  the 
solution  and  gain  faster  convergence,  while  it  can  be  seen  that  the  opposite  result  was 
achieved.  Despite  the  more  aggressive  time  step  and  expectations  of  a  more  rapid 
convergence,  the  slope  of  the  declining  residuals  when  compared  with  the  inviscid  case 
appeared  to  be  the  same. 
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Figure  11.  Mach  0.3  (Viscous)  Residual  Convergence  History 

Figures  12  and  13  illustrate  the  impact  of  the  introduction  of  viscous  terms  on 
both  the  lift  coefficient  and  the  drag  coefficient.  The  computation  of  a  more  physically 
realistic  wake,  enhanced  the  vehicle  performance  in  both  axes.  A  small,  but  notable,  drop 
in  drag  occurred  with  the  addition  of  viscosity,  which  demonstrated  that  the  addition  of 
viscous  drag  over  the  body  was  more  than  compensated  by  the  significant  reduction  in 
pressure  drag  achieved  as  the  viscous  terms  forestalled  the  onset  of  separation  at  the 
trailing  edge.  Similarly,  the  drag  coefficient,  which  continued  to  oscillate,  shifted 
upwards  to  a  new  mean  lift  coefficient  at  the  introduction  of  the  viscous  terms. 
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Figure  12.  Drag  Coefficient  vs.  Number  of  Iterations  for  Moo  =  0.3,  Viscous 
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Figure  13.  Lift  Coefficient  vs.  Number  of  Iterations  for  Mo  =  0.3,  Viscous 
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The  graphical  representation  of  the  surface  pressure  and  Mach  number  contours  in 
the  viscous  case  at  Moo  =  0.3  are  shown  in  Figure  14.  The  presence  of  Mach  number 
variation  in  the  flowfield  over  the  upper  surface  was  an  indication  of  the  boundary  layer 
development  and  propagation  into  the  flowfield.  The  lower  surface  boundary  layer  was 
visible  along  the  engine  cowl.  The  region  of  the  trailing  edge  showed  a  narrower  wake 
that  was  believed  to  be  more  representative  of  the  expected  physical  wake,  due  to  the 
presence  of  viscosity  in  the  flow.    The  computed  viscous  wake  (Figure  14)  was  smaller 
and  better  defined  than  the  wake  in  the  inviscid  case  (Figure  1 0),  which  indicated  less 
separated  flow  in  the  region  of  the  trailing. 


Figure  14.  Normalized  Surface  Pressure  and  Mach  Number  Contours,  Ho  =  0.3,  Viscous 
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3.  Subsonic  Trends 

The  convergence  histories,  force  coefficient  histories  and  FAST  images  for  all 
solutions  can  be  seen  in  Appendix  C.  Flow  solutions  for  the  remainder  of  the  subsonic 
flight  conditions  demonstrated  many  similar  characteristics.  Residual  plots  demonstrated 
the  requisite  convergence  to  a  steady  state  solution.  The  addition  of  viscous  forces  after  a 
number  of  inviscid  iterations  caused  a  spike  to  appear  in  the  residual  files  similar  to  the 
spike  seen  in  Figure  1 1 .  Beginning  at  Moo  =  0.5,  the  presence  of  the  long-period  acoustic 
resonance  dampens  out  in  the  early  iterations  leaving  only  the  shot  period  vortex- 
shedding  oscillation  evident  within  the  residual  plot.  An  additional  undamped  oscillation 
occurred  in  the  nearly  converged  viscous  solutions  at  Mo  =  0.7  (See  Figure  C-10) ,  0.8 
and  0.9.  Visual  inspection  of  the  solutions  of  these  flight  conditions  using  FAST  showed 
the  separation  and  recirculation  of  flow  in  the  wake  of  the  vehicle.  This  separation  was 
believed  to  be  the  cause  of  the  unsteadiness  in  the  residual. 

With  increasing  Mach  number  the  stagnation  point,  evidenced  by  the  Mach 
contours  seen  at  the  nose,  migrated  forward  from  a  position  aft  of  the  nose  to  a  position 
coincidence  with  the  nose  of  the  waverider  (See  Figure  C-47  to  Figure  C-60). 

Lift  and  drag  coefficient  histories  showed  very  little  difference  in  patterns  across 
the  subsonic  Mach  number  spectrum;  each  file  showed  oscillating  values  which  damped 
to  a  single  value  as  the  solution  converged.  Behaviors  of  the  lift  and  drag  coefficients, 
with  the  addition  of  viscous  terms,  were  similar  to  the  case  of  Moo  =  0.3  viscous  flow. 
Trends  in  the  values  of  these  coefficients  are  the  subject  of  the  later  section  on  Reynolds 
number  effect. 

B.         SUPERSONIC  RESULTS 

The  supersonic  solutions  had  the  advantage  of  requiring  significantly  less 
computational  time  to  obtain  a  fully  converged  solution.  The  Moo  -  4.0  solution 
converged  in  4,100  iterations,  as  opposed  to  the  160,000  iterations  required  for  Moo  =  0.5. 
In  addition  to  the  speed  at  which  the  codes  ran,  a  number  of  other  differences  between  the 
subsonic  and  supersonic  solutions  exists. 
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1.  Mach  6.0  Inviscid  Results 

Figure  13  shows  the  convergence  history  for  the  inviscid  Mo  =  6.0  solution.    This 
solution  required  some  13,000  iterations  due  to  the  small  time  step  required  to  prevent  the 
solution  from  going  unstable.  Nonetheless,  the  computational  time  was  an  order  of 
magnitude  lower  than  the  subsonic  solutions.  Note  that  both  the  long-period  oscillations, 
attributed  to  acoustic  wave  reverberation  within  the  computational  space,  and  the  short 
period  oscillations,  attributed  to  trailing  edge  vortex  shedding,  seen  in  Figures  7  and  1 1, 
were  absent  in  the  Mo  =  6.0  convergence  history. 

The  shape  of  the  convergence  history,  when  compared  to  the  subsonic  results 
(Figures  7  and  11),  was  markedly  different.  The  residuals  failed  to  drop  off  at  the  nearly 
steady  rate  seen  in  the  subsonic  convergence  histories,  but,  instead,  maintained  a  nearly 
constant  value  until  late  in  the  computation  and  finally  dropped  off  sharply.  The  reason 
for  the  convergence  history  remaining  nearly  constant  and  then  suddenly  dropping  off 
was  that,  within  the  flowfield  of  the  waverider  at  supersonic  conditions,  the  bulk  of  the 
flowfield  was  largely  unaffected,  thus  the  residuals  did  not  converge  until  the  effected 
sections  were  fully  updated.  Once  the  solution  began  to  update  the  region  close  to  the 
body,  the  solution  converged  rapidly. 
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Figure  15.  Mach  6.0  (lnviscid)  Residual  Convergence  History 

Despite  the  lack  of  a  fully  converged  solution  early  in  the  computation  the  drag 
(Figure  16)  and  lift  (Figure  17)  coefficients  rose  quickly  and  remained  nearly  steady 
throughout  the  computation.  Also  noteworthy  was  the  drop  in  the  drag  coefficient  by 
half,  to  CD  =  0.00439  from  the  Mo  =  0.3  inviscid  solution  of  CD  =  0.0194.  As  expected, 
the  lift  coefficient  at  M*,  =  6.0  is  positive.  At  Mach  6.0,  the  converged  lift  coefficient 
was  CL  =  0.018. 
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Figure  16.  Drag  Coefficient  vs.  Number  of  Iterations  for  ML  =  6.0,  Inviscid 
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Figure  17.  Lift  Coefficient  vs.  Number  of  Iterations  for  Mo  =  6.0,  Inviscid 


Figure  18  shows  the  normalized  pressure  on  the  surface  of  the  waverider  as  well 
as  the  computed  Mach  number  contours  in  the  flowfield  surrounding  the  body.  The 
pressure  differential  between  the  fore  and  aft  bodies  of  the  vehicle  appeared  less 
pronounced  than  at  Mx>  =  0.3,  however  the  scale  of  the  legend  was  changed  to 
accommodate  the  large  pressure  rise  at  the  engine  intake.  Pressures  of  nearly  14  times 
the  freestream  static  pressures  were  seen  in  this  region  of  the  flow. 


Figure  18.  Normalized  Surface  Pressure  and  Mach  Number  Contours,  Mo  =  6.0,  Inviscid 

The  Mach  contours  illustrate  the  formation  of  shocks  at  the  nose  of  the  vehicle 
and  at  the  engine  inlet.  From  the  NACA  Report  1 135  Oblique  Shock  Chart  [Ref.  13],  the 
predicted  shock  angle  for  Moo  =  6.0  flow  over  the  5.2°  wedge  of  the  vehicle  nose  is  13°. 
Although  difficult  to  measure  from  the  FAST  image,  this  appeared  to  be  the  angle  at 
which  the  shock  was  computed.  At  that  freestream  Mach  number  the  shocks  over  the 
engine  inlet  two-ramp  configuration  coalesced  into  a  single  oblique  shock.  The  engine 
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inlet  however,  when  modeled  as  a  solid  body,  dramatically  changed  the  flow  angle  and 
caused  a  distinct  shock  to  form  at  the  engine  inlet.  Here,  where  the  engine  inlet  was 
modeled  as  a  solid  body,  the  engine  inlet  shock  formed  at  a  much  higher  angle  than  the 
shocks  off  the  nose  of  the  vehicle. 


2. 


Mach  6.0  Viscous  Results 


Due  to  the  tendency  of  the  Mo  =  6.0  solution  to  go  unstable  with  minor  changes 
in  the  inputs,  the  viscous  solution  was  computed  from  zero  iterations.  This  convergence 
history  (Figure  19)  shows  that  the  viscous  solution  converged  more  slowly  than  the 
inviscid  Mao  =  6.0  solution  (Figure  15).  The  viscous  effects  required  additional  iterations 
to  allow  development  of  the  boundary  layer.  The  drag  coefficient  and  lift  coefficient 
histories  shown  in  Figures  20  and  21,  also  reflected  the  presence  of  viscous  effects  from 
the  beginning  of  the  computations,  rather  than  showing  a  change  within  their  values  when 
viscous  forces  were  added. 
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Figure  19.  Mach  6.0  (Viscous)  Residual  Convergence  History 
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Figure  20.  Drag  Coefficient  vs.  Number  of  Iterations  for  Mo  =  6.0,  Viscous 
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Figure  21.  Lift  Coefficient  vs.  Number  of  Iterations  for  Mo  =  6.0,  Viscous 
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In  Figure  22,  the  presence  of  Mach  contours  in  the  flowfield  near  the  upper 
surface  of  the  Price  waverider  was  evidence  of  the  boundary  layer  formation  on  that 
surface.  Shocks  were  more  smeared  than  in  the  inviscid  case  and  the  wake  region 
extended  over  a  larger  distance  downstream  of  the  trailing  edge. 


Figure  22.  Normalized  Surface  Pressure  and  Mach  Number  Contours,  Moo  =  6.0,  Viscous 

3.  Supersonic  Trends 

Inviscid  and  viscous  solutions  were  computed  for  MK  =  1.2  and  Moo  =  4.0,  with 
similar  results  to  the  Moo  =  6.0  solutions.  The  engine  flow-through  model  was  only 
computed  for  the  Moo  =  6.0  flight  condition.  The  results  for  additional  supersonic  cases 
are  provided  in  Appendix  C. 
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C.        ENGINE  FLOW-THROUGH  MODELING 

The  presence  of  the  extremely  high  pressure  forward  of  the  intake,  seen  in  Figures 
18  and  22,  did  not  necessarily  represent  the  true  nature  of  the  flowfield  of  an  operating 
vehicle.  To  provide  data  on  a  more  representative  configuration  engine  flow-through  was 
modeled  as  a  follow-on  test  case  to  the  M„  =  6.0  computations.  This  configuration  was 
not  modeled  for  the  subsonic  case,  although  the  impact  is  expected  to  be  less  significant 
because  the  flow  in  front  of  an  engine  in  subsonic  flight  is  expected  to  be  slowed 
significantly  at  the  diffuser  and  thus  will  present  a  surface  that  more  closely  resembles  a 
solid  body. 

The  Mach  6.0  engine  flow-through  computation  could  not  be  run  from  the 
previously  converged  M„  =  6.0  solution  because  it  not  only  required  modification  to  the 
boundary  conditions,  but  required  the  definition  of  a  new  set  of  boundaries.  How- 
through  was  extrapolated  from  the  pressure  forward  of  the  intake  boundary  and  the 
exhaust  was  modeled  as  freestream  pressure  extrapolated  flow. 

Convergence  history  for  the  ML  =  6.0  engine  flow-through  case  (Figure  23), 
shows  more  moderate  drops  in  the  residuals  than  in  the  initial  M„  =  6.0  solution  (Figure 
15).  Although  the  residuals  began  to  drop  off  between  7000  and  8000  iterations,  as  seen 
in  the  previous  Mach  6.0  solution,  the  engine  Flow-Through  model  required  nearly  twice 
as  many  iterations  to  converge,  as  the  pressure  extrapolations  from  previous  iterations 
were  computed  as  input  to  follow-on  iterations. 
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Figure  23.  Mach  6.0  (Inviscid)  Engine  Flow-Through  Convergence  History 

The  resulting  solution  was  computed  to  show  decreases  in  both  lift  and  drag 
coefficients  over  the  Mo  =  6.0  solution.  The  lift  coefficient,  calculated  in  the  inviscid 
Moo  =  6.0  solution  was  Cl  =  0.01 8,  was  found  in  the  engine  flow-through  model  to  be 
Cl  =  0.016.  The  drag  coefficient  was  found  to  be  Cd  =  0.003  for  the  engine  flow-through 
model,  whereas  previous  Mo  =  6.0  solution  (without  engine  flow-through)  computed  the 
drag  coefficient  to  be  Cd  =  0.0044.  In  this  calculation,  the  declining  lift  was  outpaced  by 
the  improvements  in  the  drag,  resulting  in  better  overall  L/D  performance  than  the 
original  calculations  demonstrated. 

As  is  illustrated  in  Figure  24,  the  pressure  at  the  intake  dropped  to  only  twice  the 
freestream  static  pressure,  as  opposed  to  the  factor  of  14  seen  in  the  initial  Mo  =  6.0 
solution.  Mach  number  contours  also  showed  that  the  resulting  shocks  occurred  closer  to 
the  body. 
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Figure  24.  Normalized  Surface  Pressure  and  Mach  Number  Contours,  Moo  =  6.0,  Engine 

Flow-Through,  Inviscid 


D. 


REYNOLDS  NUMBER  EFFECTS 


Performance  of  the  Price  waverider  across  the  spectrum  of  flight  Mach  numbers 
was  compared  using  the  lift  and  drag  coefficients  and  the  overall  performance  factor  of 
the  lift-to-drag  ratio.  The  values  for  lift  and  drag  coefficients  were  tabulated  from  the 
OVERFLOW!. 8b  output  file  containing  force  and  moment  data.  The  force  and  moment 
file  provided  the  force  and  moment  contributions  due  to  pressure  and  friction.  These 
results  are  compiled  in  Tables  2  and  3.  In  the  inviscid  case  the  contributions  due  to 
friction  were  zero  and  have  not  been  reported  in  Table  2.  In  the  viscous  cases,  the 
computed  contribution  of  the  viscous  effects  to  the  lift  and  drag  coefficients,  were,  in 
nearly  all  cases,  rounded  to  zero. 
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Mach 
Number 

CL 

CD 

L/D 

Cm 

0.3 

-0.0310 

0.0194 

-1.60 

-0.0387 

0.4 

-0.0329 

0.0171 

-1.92 

-0.0405 

0.5 

-0.0344 

0.0160 

-2.15 

-0.0424 

0.6 

-0.0377 

0.0159 

-2.37 

-0.0460 

0.7 

-0.0402 

0.0162 

-2.48 

-0.0494 

0.8 

-0.0421 

0.0178 

-2.37 

-0.0530 

0.9 

-0.0308 

0.0230 

-1.34 

-0.0463 

1.2 

0.0078 

0.0302 

0.26 

-0.0167 

2.0 

0.0312 

0.0152 

2.06 

0.0118 

4.0 

0.0224 

0.0064 

3.50 

0.0108 

6.0 

0.0184 

0.0044 

4.20 

0.0094 

Table  2.  Summary  of  Inviscid  Solution  Force  Coefficients 


Mach 
Number 


pressure 


CL 

friction 


CL  total 


CD 
friction 


Cd  total 


L/D 


Cm 


0.3 


0.0239 


0.0000 


-0.0238 


0.0169 


0.0001 


0.0170 


-1.40 


-0.0332 


0.4 


0.0203 


0.0000 


-0.0203 


0.0156 


0.0001 


0.0157         -1.29 


-0.0306 


0.5 


0.0176 


0.0000 


0.0150 


0.0001 


0.0150 


■1.17 


-0.0289 


0.6 


0.0149 


0.0000 


0.0149 


0.0001 


.0150        -0.99 


-0.0272 


0.7 


0.0109 


0.0000 


0.0152 


0.0001 


-0.71 


-0.0247 


0.8 


0.0059 


0.0000 


>.0059 


0.0163 


0.0000 


-0.36 


-0.0217 


0.9 


0.0002 


0.0000 


)IQ002 


0.0201 


0.0000 


-0.01 


-0.0186 


1.2 


-0.0108 


0.0000 


0.0108 


0.0295 


0.0000      0.029$ 


0.36 


-0.0129 


2.0 


-0.0323 


0.0000LVXK0323 


0.0155 


0.0000 


0.0,154 


2.10 


0.0126 


4.0 


-0.0240 


0.0000  mffl240 


0.0070 


0.0000 


0,0071 


3.40 


0.0115 


6.0 


-0.0214 


0.00001:;  0.0214 


0.0047 


0.0000 


).0047 


4.54 


0.0112 


Table  3.  Summary  of  Viscous  Solution  Force  Coefficients 

These  results  are  compiled  for  both  the  viscous  and  inviscid  solutions  and  shown 
graphically.  Figure  25  shows  the  lift  coefficient  as  a  function  of  Mach  number.  The 
inviscid  subsonic  lift  behavior  showed  the  detrimental  effects  of  separation,  although  it  is 
not  clearly  understood  what  specific  effect  caused  the  inviscid  lift  coefficients  to  diverge 
so  dramatically  from  the  viscous  lift  coefficients  subsonically.  In  the  supersonic  region, 
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viscous  and  inviscid  lift  coefficients  showed  very  good  agreement  with  each  other.  The 
implication  here  is  that  inviscid  solutions  closely  approximate  the  behavior  of  viscous 
flows,  and  as  such,  the  viscous  effects  may  be  negligible  in  the  initial  design  phase  of 
supersonic  vehicles.  The  viscous  solutions  consistently  computed  higher  values  for  the 
lift  coefficient,  which  demonstrated  the  beneficial  effects  of  the  attached  flow.  A 
significant  feature  of  this  configuration  was  that  the  lift  coefficient  does  not  become 
positive  until  the  flow  is  supersonic.  In  short,  the  Price  waverider,  in  its  current 
configuration,  produced  no  positive  aerodynamic  lift  until  the  freestream  Mach  number 
became  sonic.  Studies  of  the  drag  coefficient  and  pitching  moment  are  required  to 
determine  the  impact  of  not  have  positive  lift  at  zero  degrees  angle-of-attack. 

Wind  tunnel  testing  initially  performed  by  Cedrun  [Ref.  6]  provided  a  basis  of 
comparison  for  Huffs  follow-on  testing  [Ref.  7].  Huff  demonstrated  that  his  data 
duplicated  Cedrun' s  original  results.  Due  to  the  formatting  of  Cedrun' s  work,  the  lift  and 
drag  coefficients,  with  wall  and  blockage  corrections  for  the  wind  tunnel,  were  not 
provided  in  tabular  form,  so  for  greater  accuracy,  the  values  from  Huffs  work  are  shown 
in  the  comparisons  that  follow.  Since  the  data  of  both  sets  of  testing  agree,  the 
comparisons  made  here  to  Huffs  data  can  be  applied  to  Cedrun 's  work  as  well. 

The  lift  coefficient  from  Huffs  highest  speed  wind  tunnel  tests  is  plotted  on 
Figure  25  at  ML  =  0.15.  Seen  here  it  lies  slightly  below  the  extrapolated  path  of  the 
viscous  solution  of  computed  lift  coefficients,  demonstrating  that  the  wind  tunnel  test 
results  are  slightly  less  optimistic  than  the  computed  solution.  The  reason  for  this  is 
believed  to  be  the  addition  of  the  trailing  edge  cusp  on  the  computational  waverider 
model. 

The  single  engine  flow-through  configuration  lift  coefficient  is  also  plotted  on 
Figure  25  and  shows  the  slight  decrease  in  lift  at  M„o  =  6.0.  This  is  due  to  the  angle  of  the 
shock  formed  at  the  engine  inlet.  In  the  engine  flow-through  model,  this  shock  was 
slightly  more  oblique  than  the  shock  seen  in  the  solid  body  engine  model  at  M„,  =  6.0. 
The  more  oblique  shock  in  the  engine  flow-through  model  causes  a  smaller  pressure  rise 
across  the  shock,  thus  the  pressure  on  the  lower  surface  of  the  waverider  body  is  less. 
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Figure  25.  Lift  Coefficient  vs.  Mach  Number 

In  Figure  26,  the  plot  of  drag  coefficient  versus  Mach  number,  we  see  the 
transonic  drag  rise  across  M„  =  1.0  in  both  the  viscous  and  inviscid  solutions.  The 
dominant  effect  of  the  pressure  drag  over  any  viscous  forces  is  supported  by  the  close 
agreement  of  these  two  sets  of  data  through  the  entire  Mach  number  range  computed, 
0.3  <  Moo  <  6.0.  The  trough  in  drag  in  the  subsonic  region  was  an  unexpected  occurrence. 
The  subsonic  flowfield  results,  shown  in  Appendix  C,  indicate  that  the  stagnation  point 
migrated  from  a  position  slightly  aft  and  under  the  nose  of  the  vehicle  at  M„  =  0.3,  to  a 
position  at  the  nose  for  M„  =  0.5  or  0.6.  This  stagnation  point  migration  presented  a 
more  streamlined  body  to  the  freestream  flow  and  resulted  in  a  decrease  in  the  drag 
coefficient.  This  streamlining  effect  was  overcome  by  the  drag  rise  in  the  transonic 
region.  This  transonic  drag  rise  was  the  result  of  the  presence  of  regions  of  supersonic 
flow  that  had  to  shock  back  to  freestream  along  the  surface  of  the  body.  The  presence  of 
these  shocks  caused  the  early  onset  of  separation,  which  increased  the  pressure  drag  felt 
by  the  body. 
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Included  within  the  graphical  representation  of  drag  coefficients  is  the  Mo  =  0.15 
value  of  Cd  from  Huffs  wind  tunnel  tests.  This  point  lies  slightly  above  the 
extrapolation  of  the  computed  solutions.  Also  plotted  is  the  drag  coefficient  for  the  Mo  = 
6.0  engine  flow-through  configuration.  As  in  Figure  24,  the  modeling  of  engine  flow- 
through  reduced  the  shock  angle,  and  presented  a  more  streamlined  body  to  the 
freestream  flow,  which  reduced  the  drag  significantly. 
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Figure  26.  Drag  Coefficient  vs.  Mach  Number 

Overall  performance  of  the  Price  waverider  is  shown  in  Figure  27.  The  plot  of 
lift-to-drag  ratio  versus  Mach  Number  shows  the  optimum  performance  at  Moo  =  6.0. 
This  result  was  expected,  since  the  vehicle  configuration  was  optimized  for  Moo  =  6.0, 
although  this  performance  indicator  did  not  reflect  the  product  of  lift-to-drag  ratio  and 
specific  impulse.  This  study  did  not  include  the  estimation  of  engine  performance.  It 
should  also  be  noted  that,  for  purposes  of  consistency  throughout  the  study,  all  Mach 
numbers  and  Reynolds  numbers  were  computed  at  sea  level,  an  unrealistic  environment 
for  the  MU  =  6.0  flight  condition. 
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The  engine  flow-through  performance  indicated  that  the  decrease  in  drag  in  this 
configuration  outpaced  the  decline  in  lift  performance.  Also  plotted  are  the  M»  =  0.15 
lift-to-drag  ratio  for  Huffs  wind  tunnel  test  data  which  lies  in  the  plane  of  the 
extrapolated  L/D  curve  from  the  computed  solutions. 
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Figure  27.  Lift-to-Drag  Ratio  vs.  Mach  Number 

Maximum  L/D  performance  for  a  given  design  Mach  number  waverider  is  shown 
for  the  supersonic  range.  Shown  here,  Kuchemann's  maximum  on-design  performance  is 
compared  to  the  off-design  performance  data  for  the  Price  waverider.  This  estimate  was 
based  on  Kuchemann's  equation  [Ref.  3]  of  what  should  be  attainable  in  optimized 

design: 

I_3(Mro+3) 
D  M„ 

Kuchemann's  equation  is  based  on  empirical  data  and  predicts  that,  in  the  limit  as 

Mo,  goes  to  infinity,  L/D  will  asymptote  to  3.0.  At  the  design  Mach  number  of  Moo  =  6.0, 
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all  L/D  performance  values  are  shown  to  converge,  demonstrating  close  agreement 
between  the  computed  and  the  theoretical  L/D  performance  of  the  vehicle.  The  predicted 
value  of  L/D  at  Moo  =  6.0  is  4.5.  This  performance  is  closely  matched  by  the  computed 
values  of  L/D  =  4.19  for  the  inviscid  solution,  L/D  =  4.54  for  the  viscous  solution  and 
L/D  =  5.17  for  the  engine  flow-through  solution.  Note  that  both  the  viscous  solution  and 
the  engine  flow-through  solution  exceed  that  maximum  predicted  value  for  L/D. 
Anderson  [Ref.  1],  disputed  Kuchemann's  performance  limits,  estimating  that  L/Dmax 
for  supersonic  vehicles  should  be  about  50%  higher  than  Kuchemann's  predictions  for 
Mach  6.0  optimized  vehicles.  Anderson,  however,  based  that  dispute  on  a  class  of 
waverider-conflgured  vehicles  generated  using  his  viscous-optimization  routine.  While 
the  viscous  effects  have  been  taken  into  account  in  the  computations  reported  herein,  the 
Price  waverider  was  designed  using  an  inviscid  optimization  routine. 

Figure  28  shows  the  value  of  the  pitching  moment  coefficient  for  the  computed 
solutions.  Based  on  the  coordinate  system  used,  the  positive  pitching  moment  represents 
a  nose  down  pitching  force.  Without  employment  of  control  surfaces,  the  Price 
waverider  appears  to  be  statically  unstable.  The  datum  employed  here  is  the  y-axis,  the 
nose  of  the  vehicle.  The  inviscid  solution  shown  here  continues  to  demonstrate  the 
adverse  effect  of  the  non-physically  representative  wake  in  the  CFD  solution,  and  the 
subsonic  values  for  Cm  do  not  reflect  a  realistic  solution. 
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Figure  28.  Pitching  Moment  Coefficient  vs.  Mach  Number 


As  is  seen  in  Figures  25  through  28,  the  computed  force  coefficients  demonstrated 
agreement  with  the  experimental  data  from  Huffs  low-speed  wind  tunnel  testing. 
Overall,  the  modeling  of  engine  flow-through  improved  the  performance  data.  Close 
agreement  between  the  viscous  and  inviscid  solutions  in  the  supersonic  region  supported 
the  premise  that,  for  a  first  approximation,  viscous  effects  can  be  neglected  in  the  initial 
design  phases.  This  argument  cannot  be  supported  for  the  subsonic  region,  as  the  effects 
of  the  non-physically  representative  wake  cause  lift,  L/D  performance  and  pitching 
moment  parameters  to  diverge  from  the  viscous  case. 
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Computations  for  the  range  of  angles  of  attack  of -10°  <  oc  <  +10°  were 
conducted  at  MM  =  0.3,  for  viscous  conditions.  The  convergence  history,  as  well  as  lift 
and  drag  coefficient  history,  were  similar  to  those  seen  in  the  zero  degree  angle  of  attack 
solution  for  M»  =  0.3,  viscous  flow.    The  calculated  effect  on  the  lift  coefficient  is  shown 
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in  Figure  29.  Wind  tunnel  data  for  Mo  =  0.15  from  Huff's  thesis  are  shown  for 
comparison.  The  computed  values  demonstrated  a  lift  curve  slope  of  Cia  =  0.884, 
compared  to  the  wind  tunnel  test  lift  curve  slope  of  C^  =  2.526.  This  discrepancy  is 
most  likely  due  to  the  coarseness  of  the  grid.  It  is  also  possible  that  a  better  match  may 
be  obtained  by  using  a  different  turbulence  model  for  the  computations.  Cedrun's  test 
results,  although  not  plotted,  matched  HufFs  data.  Cedrun  noted  no  marked  difference  in 
the  lift  curve  slope  relative  to  Mach  number,  so  the  impact  of  the  difference  in  freestream 
Mach  number  between  the  wind  tunnel  data  and  the  computed  solutions  was  discounted 
as  the  cause  for  the  differential  in  lift  curve  slope. 
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Figure  29.  Lift  Coefficient  vs.  Angle  of  Attack 

The  onset  of  stall  was  not  predicted  to  occur  within  the  range  of  angles  of  attack 
for  which  solutions  were  computed.  Cedrun's  lowest  speed  wind  tunnel  test  (M*.  ~ 
0.003)  concluded  that  stall  occurred  at  34  degrees  AOA.  Neither  Huff  nor  Cedrun 
achieved  stall  in  higher  speed  wind  tunnel  tests,  due  to  the  total  force  vector  exceeding 
the  wind  tunnel  sting  force  constraint. 
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Figure  30  shows  the  impact  of  angle  of  attack  on  the  drag  coefficient  for  both 
wind  tunnel  data  and  computed  solutions.  The  computed  drag  coefficient  was 
consistently  lower  than  the  wind  tunnel  test  data  The  optimistic  prediction  of  drag 
coefficients  was  presumed  to  be  the  result  of  the  addition  of  the  trailing  edge  cusp.  The 
addition  of  the  cusp  also  caused  a  smaller  variation  in  drag  coefficient  over  the  range  of 
angles  of  attack,  thus  flattening  the  curve  when  compared  to  the  wind  tunnel  data. 
Minimum  drag  was  found  to  occur  at  zero  degrees  angle  of  attack. 
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Figure  30.  Drag  Coefficient  vs.  Angle  of  Attack 

A  final  comparison  was  made  between  the  wind  tunnel  data  and  the  computed 
solutions'  L/D  performance.  Shown  in  Figure  31,  the  measurements  on  the  wind  tunnel 
model  showed  a  steeper  variation  of  L/D  with  changing  angle  of  attack.  The  lift  curve 
slope  variation  was  the  largest  contributor  to  the  difference  in  L/D  performance,  however 
the  addition  of  the  trailing  edge  cusp  also  played  a  role  in  the  differences  in  performance. 
Maximum  L/D  performance  was  calculated  to  be  L/D  =  2.87,  which  occurred  at  eight 
degrees  angle  of  attack.  The  wind  tunnel  data  showed  a  maximum  of  L/D  =  3.58,  which 
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occurred  at  four  degrees  of  angle  of  attack.  This  value  may  have  continued  to  rise  had 
higher  AOA's  been  tested.  At  negative  angles  of  attack  the  computed  L/D  performance 
predicts  a  minimum  L/D  between  -A  and  -8  degrees,  closely  matching  the  wind  tunnel 
measurements,  which  measured  a  minimum  L/D  at  -6  degrees. 
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Figure  31.  Lift-to-Drag  Ratio  vs.  Angle  of  Attack 

The  normalized  surface  pressure  and  Mach  number  contours  for  the  positive  10 
degrees  angle  of  attack  case  are  shown  in  Figure  32.  The  sharp  Mach  number  contour 
gradient  showed  the  presence  of  the  boundary  layer  on  the  upper  surface  of  the  body. 
Mach  contours  outlined  the  stagnation  point  located  aft  of  the  nose  on  the  lower  surface 
and  indicated  a  region  of  separated  flow  on  the  upper  surface.  The  portion  of  the  upper 
surface  where  flow  was  separated  was  also  indicated  by  the  region  of  lower  surface 
pressure  seen  along  the  leading  edge  from  the  nose  to  mid-chord. 
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Figure  32.  Normalized  Surface  Pressure  and  Mach  Contours,  Mach  0.3  at  10°  AOA 


F. 


CFD  VALIDITY 


Based  on  the  data  presented  above,  the  CFD  appeared  to  calculate  solutions, 
which  represented  realistic  flowfields.  While  computational  limitations  prohibited 
calculation  of  the  flowfield  solution  below  Mo  =  0.3,  the  computed  lift  and  drag 
coefficients  and  L/D  for  0.3  <  Mo  <  6.0  trended  toward  the  values  obtained  in  the  Moo  = 
0.15  tests  conducted  by  Huff  and  Cedrun.  While  precise  comparison  was  not  made,  due 
to  the  mismatch  of  Mach  numbers  in  the  CFD  and  wind  tunnel  test,  it  appeared  that  the 
CFD  trends  were  slightly  optimistic.  The  optimism  of  the  CFD  calculations  for  lift  and 
drag  was  believed  to  be  the  result  of  the  addition  of  the  trailing  edge  cusp  to  the  CFD 
model. 

Angle  of  attack  sweeps  showed  disagreement  in  all  computed  performance 
parameters.  While  some  of  the  discrepancy  was  attributable  to  the  addition  of  the  trailing 
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edge  cusp,  the  significant  difference  in  the  lift  curve  slope  is  believed  to  be  the  result  of 
the  grid  coarseness  and  turbulence  model.  Despite  disagreement  at  both  positive  and 
negative  angles  of  attack,  the  solutions  at  zero  degrees  angle  of  attack  were  in  close 
agreement.  For  that  reason,  the  zero  degree  AOA  data  were  still  considered  valid. 

The  flight  conditions  considered  herein  represent  the  Mach  number  and  Reynolds 
number  calculated  for  the  15-inch  model  at  sea  level.  To  better  approximate  the 
boundary  layer  development,  Reynolds  numbers  should  be  chosen  to  reflect  the  length  of 
the  Price  waverider,  rather  than  the  model  root  mean  chord. 


G.         CFD  LIMITATIONS 

For  purposes  of  this  study,  a  broad  range  of  Mach  numbers  was  examined, 
however  the  extent  to  which  existing  data  could  be  matched  was  limited  by  the  flow 
solver.  Flow  solutions  were  not  obtained  for  Mach  numbers  below  Moo  =  0.3,  whereas 
wind  tunnel  test  data  was  not  available  above  Moo  =  0.15.  This  mismatch  prevented  the 
validation  of  CFD  results  with  corresponding  wind  tunnel  data. 

Difficulty  in  creating  a  grid  that  correctly  modeled  the  blunt  trailing  edge,  led  to 
the  addition  of  the  trailing  edge  cusp  to  the  CFD  model.  This  difference  in  configuration 
is  believed  to  have  caused  the  lift  and  drag  coefficients  to  be  more  optimistic  than  the 
wind  tunnel  test  data. 

The  use  of  the  elliptic  solver,  GRIDGEN3D,  to  resolve  the  negative  cell  volumes 
in  the  grid  also  tended  to  pull  grid  points  away  from  the  body  of  the  waverider.  Due  the 
migration  of  grid  points  away  from  the  body  of  the  waverider,  the  boundary  layer 
calculations,  while  seen  in  the  Mach  number  contours,  may  not  be  representative  of 
actual  boundary  layer  size. 
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V.         CONCLUSIONS 

The  computational  fluid  dynamic  analysis  of  the  Price  waverider  was 
demonstrated  to  produce  valid  results  for  subsonic  and  supersonic  flowfield 
computations.  The  prediction  of  lift  and  drag  coefficients  and  L/D  closely  resembled  the 
data  obtained  in  wind  tunnel  testing.  Although  the  lift  and  drag  coefficients  appeared 
slightly  more  optimistic  than  the  wind  tunnels  tests,  this  was  likely  the  result  of  the 
addition  of  the  cusped  trailing  edge.  While  that  modification  was  made  to  facilitate  the 
grid  generation,  the  physical  impact  was  to  reduce  the  pressure  drag  and  increase  the  lift 
through  the  expansion  over  the  upper  surface  trailing  edge.  The  differences  in  the 
computed  and  measured  values  were  in  the  same  order  of  magnitude  as  the  tunnel  wall 
corrections  applied  to  the  wind  tunnel  results. 

The  use  of  computational  methods  provided  the  opportunity  to  expand  the  body  of 
knowledge  available  defining  the  behavior  of  flow  surrounding  the  Price  waverider  over 
a  broad  range  of  flight  conditions. 

The  angle  of  attack  sweep  calculation,  while  showing  similar  trends  in 
performance  parameters,  did  not  achieve  very  good  agreement,  except  at  zero  degrees 
angle  of  attack.  The  lift  curve  slope  showed  the  largest  discrepancy  when  compared  to 
the  wind  tunnel  data.  This  was  believed  to  be  the  result  of  the  grid  coarseness.  Further 
study  is  required  to  determine  whether  grid  refinements  or  use  of  a  different  turbulence 
model  would  improve  agreement. 

Instability  of  the  solution  at  high  Mach  numbers  precluded  the  extension  of  CFD 
analysis  to  flight  conditions  above  the  design  condition.  Time  steps  required  to  maintain 
a  stable  solution  at  M„  =  6.0  were  an  order  of  magnitude  smaller  than  any  other 
computation.  Similarly,  the  excessive  computational  time  required  to  compute  the  low 
subsonic  solutions,  prohibits  the  use  of  OVERFLOW  1.8b  to  speeds  lower  than  M*,  =  0.3. 

The  computed  flow  solutions  demonstrated  several  areas  of  concern.  The  first 
concern  was  the  lack  of  positive  lift  generation  in  subsonic  flight  at  zero  degrees  angle  of 
attack.  With  a  tendency  to  pitch  up  at  this  angle  of  attack  present,  the  risk  of  flying  at  a 
slightly  nose  up  attitude,  to  gain  positive  lift,  is  that  this  attitude  will  exacerbate  the 
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pitching  moment.  The  second  area  of  concern  was  the  transonic  drag  rise  and  the  effect 
of  the  drag  trough  in  the  subsonic  regime  on  the  specific  excess  power  curves. 

A  number  of  opportunities  exist  to  improve  the  body  of  knowledge  regarding  the 
Price  waverider  and  in  the  extension  of  the  existing  analysis. 

A.        GRID  REFINEMENTS 

The  existing  grid  provided  a  computational  space  that  allowed  the  computation  of 
reasonable  solutions,  however,  there  were  a  number  of  areas  that  needed  improvement. 
The  existing  grid,  after  using  the  elliptical  solver  to  resolve  the  presence  of  negative  cell 
volumes,  failed  to  maintain  the  initially  specified  grid  spacing  near  the  body  except  at  the 
centerline  where  all  points  were  held  as  fixed  points.  This  spacing  contributed  to 
inaccuracies  in  the  computation  of  viscous  drag  insofar  as  the  velocity  gradient  in  the 
boundary  layer  was  dependent  on  the  existence  of  sufficient  points  near  the  body  to 
capture  the  gradient  correctly.  Additionally,  the  lift  curve  slope  in  the  alpha  sweep  did 
not  mirror  the  wind  tunnel  results.  Grid  density  was  believed  to  be  the  reason  for  this 
disagreement. 

The  addition  of  the  cusp  at  the  trailing  edge  of  the  Price  waverider,  while  a  likely 
design  change  in  the  real  world,  did  not  reflect  the  vehicle  as-tested.  It  was  seen  that  the 
lift  and  drag  forces  computed  over  the  existing  grid  were  optimistic  when  compared  to 
the  wind  tunnel  test  data.  Use  of  multi-block  grid  generation  techniques  may  improve  the 
opportunity  to  match  the  grid  geometry  to  the  wind  tunnel  model 

It  should  be  noted  that  the  leading  edge  of  the  computational  model  was 
represented  by  a  single  point  line,  with  no  leading  edge  radius.  This  contributed  to  a 
departure  from  both  the  model,  as-tested,  and  the  real-world  vehicle  requirements.  While 
the  leading  edge  radius  of  the  Price  waverider  was  not  included  in  the  design  criteria,  it 
can  be  safely  assumed  that  the  leading  edge  would  have  to  have  some  bluntness.  An 
infinitely  thin  leading  is  unrealistic  to  build,  and  would  not  be  able  to  withstand  the  high 
heating  conditions  produced  at  hypersonic  speeds.  This  lack  of  bluntness  at  the  leading 
edge  contributed  to  unrealistic  boundary  layer  development. 
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In  general,  the  grid  was  fairly  coarse,  as  only  45  points  along  the  length  of  the 
waverider  body  were  used  during  the  computations.  Grid  dependence  could  be  seen  in 
the  unevenness  of  the  Mach  contours  in  several  of  the  solutions. 

B.         OPPORTUNITIES  FOR  FURTHER  TESTING 

The  use  of  a  multi-grid  computational  model  would  enable  the  closer  examination 
of  some  aspects  of  the  flowfield.  This  method  could  also  be  used  to  make  modifications 
to  the  basic  waverider,  to  include  cockpit,  vertical  tails  or  high-lift  devices  to  gain  a  better 
understanding  of  the  aerodynamic  impacts  of  varying  configurations.  Included  in  this  list 
is  the  modeling  of  the  engine  as  a  thrust  producing  element  rather  than  as  a  solid  body  or 
a  flow-through  duct.  As  with  the  employment  of  a  cusped  trailing  edge,  a  number  of 
configuration  changes  or  proposed  changes  will  require  analysis  before  the  feasibility  of 
the  Price  waverider,  or  a  derivative,  can  be  fully  assessed. 

This  study  did  not  examine  the  effects  of  yaw  angles  or  the  combinations  of 
angle-of-attack  and  yaw  on  the  flowfield.  The  opportunity,  therefore,  exists  to  further  the 
understanding  of  the  flight  characteristics  of  the  Price  waverider  at  those  off-design  flight 
conditions. 

To  date,  supersonic  wind  tunnel  testing  has  not  been  accomplished  for  the  Price 
waverider.  A  study  that  combined  supersonic  wind  tunnel  test  data,  with  the  results 
provided  herein  would  serve  to  validate  the  CFD  code  further. 
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APPENDIX  A.  WAVERIDER  SURFACE  GRID  POINTS 
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0.07329 

0.76178 

0.41431 

0.07592 

34 

0.81896 

0.00000 

0.10925 

0.81896 

0.01870 

0.10925 

0.81896 

0.03735 

0.10925 

0.81896 

0.05595 

0.10925 

0.81896 

0.07450 

0.10925 

0.81896 

0.07801 

0.10925 

0.81896 

0.07818 

0.10925 

0.81896 

0.07828 

0.07884 

0.81896 

0.08959 

0.07862 

0.81896 

0.10727 

0.07815 

0.81896 

0.12483 

0.07762 

0.81896 

0.14225 

0.07694 

0.81896 

0.15952 

0.07620 

0.81896 

0.17663 

0.07542 

0.81896 

0.19356 

0.07463 

0.81896 

0.21029 

0.07382 

0.81896 

0.22681 

0.07305 

0.81896 

0.24310 

0.07235 

0.81896 

0.25914 

0.07169 

0.81896 

0.27494 

0.07119 

0.81896 

0.28884 

0.07087 

0.81896 

0.29049 

0.07088 

0.81896 

0.30574 

0.07073 

0.81896 

0.32070 

0.07079 

0.81896 

0.33535 

0.07111 

0.81896 

0.34972 

0.07177 

0.81896 

0.35778 

0.07224 

0.81896 

0.36374 

0.07270 

0.81896 

0.37744 

0.07400 

0.81896 

0.39081 

0.07568 

0.81896 

0.39248 

0.07591 

0.81896 

0.40385 

0.07778 

0.81896 

0.41655 

0.08030 

0.81896 

0.42889 

0.08323 

35 

0.84998 

0.00000 

0.10925 

0.84998 

0.01870 

0.10925 

0.84998 

0.03735 

0.10925 

68 


0.84998 

0.04887 

0.00262 

0.84998 

0.06515 

0.00362 

0.84998 

0.06824 

0.00382 

0.84998 

0.06839 

0.00383 

0.84998 

0.06854 

0.00384 

0.84998 

0.08143 

0.00470 

0.84998 

0.09770 

0.00582 

0.84998 

0.11397 

0.00711 

0.84998 

0.13022 

0.00848 

0.84998 

0.14647 

0.01000 

0.84998 

0.16269 

0.01165 

0.84998 

0.17890 

0.01350 

0.84998 

0.19509 

0.01553 

0.84998 

0.21125 

0.01777 

0.84998 

0.22737 

0.02024 

0.84998 

0.24345 

0.02297 

0.84998 

0.25948 

0.02598 

0.84998 

0.27377 

0.02893 

0.84998 

0.27546 

0.02930 

0.84998 

0.29135 

0.03291 

0.84998 

0.30717 

0.03690 

0.84998 

0.32290 

0.04125 

0.84998 

0.33851 

0.04599 

0.84998 

0.34740 

0.04888 

0.84998 

0.35399 

0.05110 

0.84998 

0.36934 

0.05667 

0.84998 

0.38451 

0.06264 

0.84998 

0.38643 

0.06344 

0.84998 

0.39950 

0.06904 

0.84998 

0.41431 

0.07593 

0.84998 

0.42889 

0.08323 

0.84998 

0.43654 

0.08728 

36 

0.87780 

0.00000 

0.00001 

0.87780 

0.01629 

0.00083 

0.87780 

0.03258 

0.00169 

0.87780 

0.04887 

0.00263 

0.87780 

0.06515 

0.00363 

0.87780 

0.06824 

0.00382 

0.87780 

0.06839 

0.00383 

0.87780 

0.06854 

0.00384 

0.87780 

0.08143 

0.00468 

0.87780 

0.09770 

0.00585 

0.87780 

0.11397 

0.00711 

0.87780 

0.13022 

0.00845 

0.87780 

0.14647 

0.00999 

0.87780 

0.16270 

0.01166 

0.87780 

0.17890 

0.01348 

0.84998 

0.05595 

0.10925 

0.84998 

0.07450 

0.10925 

0.84998 

0.07801 

0.10925 

0.84998 

0.07818 

0.10925 

0.84998 

0.07828 

0.08202 

0.84998 

0.08995 

0.08180 

0.84998 

0.10769 

0.08133 

0.84998 

0.12530 

0.08070 

0.84998 

0.14279 

0.08001 

0.84998 

0.16013 

0.07928 

0.84998 

0.17729 

0.07846 

0.84998 

0.19429 

0.07767 

0.84998 

0.21106 

0.07680 

0.84998 

0.22764 

0.07600 

0.84998 

0.24398 

0.07527 

0.84998 

0.26008 

0.07460 

0.84998 

0.27592 

0.07406 

0.84998 

0.28986 

0.07371 

0.84998 

0.29150 

0.07368 

0.84998 

0.30679 

0.07348 

0.84998 

0.32179 

0.07352 

0.84998 

0.33648 

0.07381 

0.84998 

0.35085 

0.07439 

0.84998 

0.35895 

0.07486 

0.84998 

0.36491 

0.07529 

0.84998 

0.37863 

0.07654 

0.84998 

0.39202 

0.07817 

0.84998 

0.39369 

0.07840 

0.84998 

0.40506 

0.08020 

0.84998 

0.41776 

0.08266 

0.84998 

0.43012 

0.08556 

0.84998 

0.43654 

0.08728 

36 

0.87780 

0.00000 

0.10925 

0.87780 

0.01870 

0.10925 

0.87780 

0.03735 

0.10925 

0.87780 

0.05595 

0.10925 

0.87780 

0.07450 

0.10925 

0.87780 

0.07801 

0.10925 

0.87780 

0.07818 

0.10925 

0.87780 

0.07828 

0.07371 

0.87780 

0.08906 

0.07377 

0.87780 

0.10668 

0.07368 

0.87780 

0.12420 

0.07350 

0.87780 

0.14161 

0.07327 

0.87780 

0.15887 

0.07292 

0.87780 

0.17600 

0.07253 

0.87780 

0.19298 

0.07221 
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0.87780 

0.19509 

0.01553 

0.87780 

0.21125 

0.01778 

0.87780 

0.22737 

0.02025 

0.87780 

0.24346 

0.02298 

0.87780 

0.25949 

0.02599 

0.87780 

0.27377 

0.02894 

0.87780 

0.27546 

0.02931 

0.87780 

0.29136 

0.03294 

0.87780 

0.30717 

0.03691 

0.87780 

0.32290 

0.04126 

0.87780 

0.33851 

0.04598 

0.87780 

0.34740 

0.04889 

0.87780 

0.35399 

0.05112 

0.87780 

0.36934 

0.05668 

0.87780 

0.38451 

0.06265 

0.87780 

0.38643 

0.06344 

0.87780 

0.39950 

0.06904 

0.87780 

0.41431 

0.07593 

0.87780 

0.42890 

0.08323 

0.87780 

0.43654 

0.08729 

0.87780 

0.44325 

0.09098 

36 

0.93164 

0.00000 

0.00002 

0.93164 

0.01629 

0.00084 

0.93164 

0.03258 

0.00171 

0.93164 

0.04887 

0.00263 

0.93164 

0.06515 

0.00363 

0.93164 

0.06824 

0.00383 

0.93164 

0.06839 

0.00384 

0.93164 

0.06854 

0.00385 

0.93164 

0.08143 

0.00469 

0.93164 

0.09770 

0.00583 

0.93164 

0.11397 

0.00711 

0.93164 

0.13022 

0.00847 

0.93164 

0.14647 

0.01000 

0.93164 

0.16269 

0.01165 

0.93164 

0.17890 

0.01350 

0.93164 

0.19509 

0.01553 

0.93164 

0.21125 

0.01777 

0.93164 

0.22737 

0.02024 

0.93164 

0.24345 

0.02298 

0.93164 

0.25949 

0.02599 

0.93164 

0.27377 

0.02893 

0.93164 

0.27546 

0.02930 

0.93164 

0.29136 

0.03293 

0.93164 

0.30718 

0.03691 

0.93164 

0.32291 

0.04127 

0.93164 

0.33852 

0.04600 

0.87780 

0.20974 

0.07173 

0.87780 

0.22633 

0.07134 

0.87780 

0.24271 

0.07103 

0.87780 

0.25886 

0.07081 

0.87780 

0.27478 

0.07070 

0.87780 

0.28880 

0.07076 

0.87780 

0.29044 

0.07075 

0.87780 

0.30584 

0.07098 

0.87780 

0.32095 

0.07142 

0.87780 

0.33573 

0.07202 

0.87780 

0.35028 

0.07308 

0.87780 

0.35844 

0.07373 

0.87780 

0.36448 

0.07433 

0.87780 

0.37837 

0.07599 

0.87780 

0.39188 

0.07789 

0.87780 

0.39356 

0.07814 

0.87780 

0.40606 

0.08221 

0.87780 

0.41863 

0.08436 

0.87780 

0.43089 

0.08703 

0.87780 

0.43728 

0.08867 

0.87780 

0.44325 

0.09098 

36 

0.93164 

0.00000 

0.10925 

0.93164 

0.01870 

0.10925 

0.93164 

0.03735 

0.10925 

0.93164 

0.05595 

0.10925 

0.93164 

0.07450 

0.10925 

0.93164 

0.07801 

0.10925 

0.93164 

0.07818 

0.10925 

0.93164 

0.07828 

0.05684 

0.93164 

0.08723 

0.05716 

0.93164 

0.10457 

0.05778 

0.93164 

0.12184 

0.05822 

0.93164 

0.13906 

0.05875 

0.93164 

0.15616 

0.05915 

0.93164 

0.17319 

0.05970 

0.93164 

0.19008 

0.06011 

0.93164 

0.20685 

0.06063 

0.93164 

0.22350 

0.06129 

0.93164 

0.23998 

0.06199 

0.93164 

0.25628 

0.06281 

0.93164 

0.27242 

0.06382 

0.93164 

0.28666 

0.06482 

0.93164 

0.28835 

0.06496 

0.93164 

0.30406 

0.06632 

0.93164 

0.31957 

0.06796 

0.93164 

0.33484 

0.06987 

0.93164 

0.34986 

0.07209 
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0.93164 

0.34740 

0.04887 

0.93164 

0.35400 

0.05113 

0.93164 

0.36933 

0.05666 

0.93164 

0.38451 

0.06265 

0.93164 

0.38642 

0.06341 

0.93164 

0.39951 

0.06906 

0.93164 

0.41431 

0.07593 

0.93164 

0.42890 

0.08323 

0.93164 

0.43654 

0.08729 

0.93164 

0.44325 

0.09098 

0.93816 

0.45582 

0.09828 

38 

1.00000 

0.00000 

0.00001 

1.00000 

0.01629 

0.00083 

1.00000 

0.03258 

0.00170 

1.00000 

0.04887 

0.00263 

1.00000 

0.06515 

0.00362 

1.00000 

0.06824 

0.00381 

1.00000 

0.06839 

0.00382 

1.00000 

0.06854 

0.00383 

1.00000 

0.08143 

0.00468 

1.00000 

0.09770 

0.00584 

1.00000 

0.11397 

0.00710 

1.00000 

0.13022 

0.00847 

1.00000 

0.14647 

0.00999 

1.00000 

0.16269 

0.01165 

1.00000 

0.17890 

0.01349 

1.00000 

0.19509 

0.01552 

1.00000 

0.21125 

0.01777 

1.00000 

0.22737 

0.02024 

1.00000 

0.24345 

0.02298 

1.00000 

0.25949 

0.02599 

1.00000 

0.27377 

0.02894 

1.00000 

0.27546 

0.02930 

1.00000 

0.29136 

0.03294 

1.00000 

0.30718 

0.03692 

1.00000 

0.32290 

0.04127 

1.00000 

0.33851 

0.04600 

1.00000 

0.34740 

0.04889 

1.00000 

0.35400 

0.05113 

1.00000 

0.36934 

0.05668 

1.00000 

0.38452 

0.06265 

1.00000 

0.38643 

0.06344 

1.00000 

0.39951 

0.06907 

1.00000 

0.41431 

0.07593 

1.00000 

0.42890 

0.08323 

1 .00000 

0.43654 

0.08729 

1.00000 

0.44325 

0.09099 

0.93164 

0.35835 

0.07352 

0.93164 

0.36458 

0.07456 

0.93164 

0.37909 

0.07753 

0.93164 

0.39323 

0.08068 

0.93164 

0.39505 

0.08121 

0.93164 

0.40798 

0.08607 

0.93164 

0.42031 

0.08764 

0.93164 

0.43239 

0.08988 

0.93164 

0.43870 

0.09135 

0.93164 

0.44446 

0.09322 

0.93164 

0.45582 

0.09828 

38 

1.00000 

0.00000 

0.03000 

1.00000 

0.01694 

0.03000 

1.00000 

0.03384 

0.03000 

1.00000 

0.05069 

0.03000 

1.00000 

0.06749 

0.03000 

1.00000 

0.07067 

0.03000 

1.00000 

0.07116 

0.03371 

1.00000 

0.07132 

0.03372 

1.00000 

0.08473 

0.03452 

1.00000 

0.10164 

0.03560 

1.00000 

0.11854 

0.03676 

1.00000 

0.13542 

0.03803 

1.00000 

0.15227 

0.03943 

1.00000 

0.16910 

0.04097 

1.00000 

0.18590 

0.04268 

1.00000 

0.20266 

0.04457 

1.00000 

0.21938 

0.04665 

1.00000 

0.23605 

0.04897 

1.00000 

0.25265 

0.05154 

1.00000 

0.26920 

0.05438 

1.00000 

0.28392 

0.05717 

1.00000 

0.28566 

0.05752 

1.00000 

0.30203 

0.06098 

1.00000 

0.31830 

0.06479 

1.00000 

0.33445 

0.06895 

1.00000 

0.35047 

0.07351 

1.00000 

0.35959 

0.07630 

'l  .00000 

0.36634 

0.07847 

1.00000 

0.38205 

0.08384 

1.00000 

0.39757 

.  0.08966 

1.00000 

0.39952 

0.09042 

1.00000 

0.41045 

0.09102 

1.00000 

0.42245 

0.09180 

1.00000 

0.43429 

0.09349 

1.00000 

0.44052 

0.09475 

1.00000 

0.44599 

0.09608 
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1.00000 

0.45736 

0.09918 

1.00000 

0.47120 

0.10780 

14 

1.01500 

0.00000 

0.01500 

1.01500 

0.01694 

0.01541 

1.01500 

0.03258 

0.01585 

1.01500 

0.07000 

0.01692 

1.01500 

0.10000 

0.02070 

1.01500 

0.13250 

0.02330 

1.01500 

0.18500 

0.02859 

1.01500 

0.25500 

0.03876 

1.01500 

0.30400 

0.04895 

1.01500 

0.32290 

0.05407 

1.01500 

0.38300 

0.07324 

1.01500 

0.42600 

0.08751 

1.01500 

0.44450 

0.09353 

1.01500 

0.47120 

0.10780 

1.00000 

0.45756 

0.09954 

1.00000 

0.47120 

0.10780 

14 
1.01500 

0.00000 

0.01500 

1.01500 

0.01694 

0.01541 

1.01500 

0.03258 

0.01585 

1.01500 

0.07000 

0.01692 

1.01500 

0.10000 

0.02070 

1.01500 

0.13250 

0.02330 

1.01500 

0.18500 

0.02859 

1.01500 

0.25500 

0.03876 

1.01500 

0.30400 

0.04895 

1.01500 

0.32290 

0.05407 

1.01500 

0.38300 

0.07324 

1.01500 

0.42600 

0.08751 

1.01500 

0.44450 

0.09353 

1.01500 

0.47120 

0.10780 
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APPENDIX  B.  SAMPLE  INPUT  FILES 


Mach  3.0  (Inviscid)  Input  file,  Start 


CHIMRA= 

.F.,    NSTEPS=200 

RESTRT=  .F.,    NSAVE 

NQT 

o, 

$END 

$FLOINP 

ALPHA  = 

0,   FSMACH  =  0.3, 

REY    =  3.19E6,   TINF 

520.000, 

$END 

$VARGAM 

$END 

$GRDNAM 

NAME  = 

'Waverider  Half-Body,  42x45x50  grid', 

$END 

$NITERS 

$END 

$METPRM 

IRHS   = 

0,    ILHS   =  2, 

IDISS   =  2, 

$END 

$TIMACU 

DT 

0.01,  ITIME=  1, 

TFOSO  =  1.00,  CFLMIN=0.0, 

$END 

$SMOACU 

ISPECJ= 

2,   DIS2J  =  0.0, 

DIS4J  =  0.2, 

ISPECK= 

2,   DIS2K  =  0.0, 

DIS4K  =  0.2, 

ISPECL= 

2,   DIS2L  =  0.0, 

DIS4L  =  0.2, 

SMOO   = 

1.00, 

EPSE   = 

0.35, 

$END 

$VISINP 

VISCJ  = 

.F. ,   VISCK  =  .F. 

,   VISCL  =  .F. , 

NTURB  = 

o, 

ITTYP  = 

1/ 

ITDIR  = 

3, 

JTLS   = 

1, 

JTLE   - 

42, 

KTLS   = 

If 

t 

KTLE   = 

45, 

LTLS   = 

1, 

LTLE   = 

50, 

TLPAR1= 

0.3, 

$END 

=10, 
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$BCINP 

NBC 

= 

6, 

IBTYP 

= 

1, 

14, 

12, 

54, 

12, 

32, 

IBDIR 

= 

3, 

2, 

-If 

-2, 

If 

-3, 

JBCS 

= 

1, 

1, 

42, 

1, 

1, 

If 

JBCE 

= 

42, 

42, 

42, 

42, 

1, 

42, 

KBCS 

= 

1, 

1, 

1, 

45, 

1, 

1, 

KBCE 

= 

45, 

1, 

45, 

45, 

45, 

45, 

LBCS 

= 

1, 

1, 

If 

1, 

1, 

50, 

LBCE 

= 

1, 

49, 

50, 

50, 

50, 

50, 

$END 

$SCEINP 

. 

$END 
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M ach  3.0  (Viscous)  Input  File,  Restart 

$ GLOBAL 

CHIMRA=  .F. ,    NSTEPS=10000,    RESTRT=  .  T . ,    NSAVE 
=100,  NQT    =  0, 

$END 
$FLOINP 

ALPHA  =  0,   FSMACH  =0.3,   REY    =  3.19E6,   TINF   = 
520.000, 

$END 
$VARGAM 

$END 
$GRDNAM 

NAME  =  'Waverider  Half-Body,  4.2x45x50  grid', 
$END 
$NITERS 

$END 
$METPRM 

IRHS   =  0,    ILHS   =  2,    IDISS   =  2, 

$END 
$TIMACU 

DT     =0.07,  ITIME=  1,   TFOSO  =  1.00,  CFLMIN=0.01, 

$END 
$SMOACU 

ISPECJ=  2,   DIS2J  =  0.0,   DIS4J  =  0.2, 

ISPECK=  2,   DIS2K  =  0.0,   DIS4K  =  0.2, 

ISPECL=  2,   DIS2L  =  0.0,   DIS4L  =  0.2, 

SMOO   =  1.00, 


VISCL 


EPSE 

= 

0.35, 

$END 

$VISINP 

VISCJ 

= 

.T. ,   VISCK  = 

NTURB 

= 

1, 

ITTYP 

= 

1, 

ITDIR 

= 

3, 

JTLS 

= 

1, 

JTLE 

= 

42, 

KTLS 

= 

1, 

KTLE 

= 

45, 

LTLS 

= 

1, 

LTLE 

= 

50, 

TLPAR] 

0.3, 

$END 

$BCINP 

NBC 

= 

6, 
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IBTYP 

= 

5, 

14, 

12, 

54, 

12, 

32, 

IBDIR 

= 

3, 

2, 

-1, 

-2, 

1, 

-3, 

JBCS 

= 

1, 

1, 

42, 

1. 

1. 

1, 

JBCE 

= 

42, 

42, 

42, 

42, 

1, 

42, 

KBCS 

= 

1, 

1, 

1, 

45, 

1, 

1, 

KBCE 

= 

45, 

1, 

45, 

45, 

45, 

45, 

LBCS 

= 

1, 

1, 

1, 

1, 

1, 

50, 

LBCE 

= 

1, 

49, 

50, 

50, 

50, 

50, 

$END 

$SCEINP 

$END 
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Mach  6.0  (Inviscid)  Input  file,  Start 

$GLOBAL 

CHIMRA=  .F.f    NSTEPS=100,    RESTRT=  .F.,    NSAVE  =10, 

NQT    =  0, 

$END 
$FLOINP 

ALPHA  =  0,   FSMACH  =  6.0,   REY    =  6.39E7,   TINF   = 
520.000, 

$END 
$VARGAM 

$END 
$GRDNAM 

NAME  =  'Waverider  Half-Body,  42x45x50  grid', 
$END 
$NITERS 

$END 
$METPRM 

IRHS   =  0,    ILHS   =  2,    IDISS   =  2, 

$END 
$TIMACU 

DT     =  0.001,  ITIME=  1,   TFOSO  =  1.00,  CFLMIN=0.0, 

$END 
$SMOACU 

ISPECJ=  2,   DIS2J  =  2.0,   DIS4J  =  0.2, 

ISPECK=  2,   DIS2K  =  2.0,   DIS4K  =  0.2, 

ISPECL=  2,   DIS2L  =  2.0,   DIS4L  =  0.2, 

SMOO   =  1.00, 

EPSE   =  0.35, 

$END 
$VISINP 

VISCJ  =  .F.,   VISCK  =  .F.,   VISCL  =  .F., 

NTURB  =    0, 

ITTYP  =    1, 

ITDIR  =    3, 

JTLS   =    1, 

JTLE   -   42, 

KTLS   =    1, 

KTLE   -   45, 

LTLS   =    1, 

LTLE   =50, 

TLPAR1=  0.3, 

$END 
$BCINP 

NBC    =    6, 

IBTYP  =      1,  14,    12,    54,    12,   32, 
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IBDIR 

= 

3, 

2, 

-1, 

-2, 

1, 

-3, 

JBCS 

= 

1, 

1, 

42, 

If 

1, 

1, 

JBCE 

= 

42, 

42, 

42, 

42, 

1, 

42, 

KBCS 

= 

1, 

1, 

1/ 

45, 

1, 

1, 

KBCE 

= 

45, 

1, 

45, 

45, 

45, 

45, 

LBCS 

= 

1, 

1, 

1, 

1, 

1, 

50, 

LBCE 

= 

1, 

50, 

50, 

50, 

50, 

50, 

$END 

$SCEINP 

$END 
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Mach  6.0  (Viscous)  File  input,  Restart 

$ GLOBAL 

CHIMRA=  .F. ,    NSTEPS=4000,    RESTRT=  .T.,    NSAVE 
=  100, 

NQT    =  0, 

$END 
$FLOINP 

ALPHA  =  0,   FSMACH  =  6.0,   REY    =  6.39E7,   TINF   = 
520.000, 

$END 
$VARGAM 

$END 
$GRDNAM 

NAME  =  'Waverider  Half-Body,  42x45x50  grid', 
$END 
$NITERS 

$END 
$METPRM 

IRHS   =  0,    ILHS   =  2,    IDISS   =  2, 

$END 
$TIMACU 

DT     =  0.005,  ITIME=  1,   TFOSO  =  1.00,  CFLMIN=0.01, 

$END 
$SMOACU 

ISPECJ=  2,   DIS2J  =  2.0,   DIS4J  =  0.2, 

ISPECK=  2,   DIS2K  =  2.0,   DIS4K  =  0.2, 

ISPECL=  2,   DIS2L  =  2.0,   DIS4L  =  0.2, 

SMOO   =  1.00, 


VISCK  =  .F.,   VISCL  = 


EPSE 

= 

0.35, 

$END 

$VISINP 

VISCJ 

= 

.T., 

NTURB 

= 

1, 

ITTYP 

= 

1, 

ITDIR 

= 

3, 

JTLS 

= 

1, 

JTLE 

= 

42, 

KTLS 

= 

1, 

KTLE 

= 

45, 

LTLS 

= 

1, 

LTLE 

= 

50, 

TLPAR] 

0.3, 

$END 
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$BCINP 
NBC 
IBTYP 
IBDIR 
JBCS 
JBCE 
KBCS 
KBCE 
LBCS 
LBCE 
$END 

$SCEINP 

$END 


6, 

5, 

14, 

12, 

54, 

12, 

32, 

3, 

2, 

-1, 

-2, 

1, 

-3, 

1, 

1, 

42, 

1, 

1, 

1, 

42, 

42, 

42, 

42, 

1, 

42, 

1, 

1, 

1, 

45, 

1, 

1, 

45, 

1, 

45, 

45, 

45, 

45, 

1, 

1, 

1, 

1, 

1, 

50, 

1, 

50, 

50, 

50, 

50, 

50, 
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M ach  6.0  (Inviscid)  Engine  Flow-Through  model  Input  file,  Start 

$GLOBAL 

CHIMRA=  .F.f    NSTEPS=100,    RESTRT=  .F.,    NSAVE  =10, 

NQT    =  0, 

$END 
$FLOINP 

ALPHA  =  0,   FSMACH  =  6.0,   REY    =  6.3  9E7,   TINF   = 
520.000, 

$END 
$VARGAM 

$END 
$GRDNAM 

NAME  =  'Waverider  Half-Body,  42x45x50  grid', 
$END 
$NITERS 

$END 
$METPRM 

IRHS   =  0,    ILHS   =  2,    IDISS   =  2, 

$END 
$TIMACU 

DT     =  0.001,  ITIME=  1,   TFOSO  =  1.00,  CFLMIN=0.0, 

$END 
$SMOACU 

ISPEGJ=  2,   DIS2J  =  2.0,   DIS4J  =  0.2, 

ISPECK=  2,   DIS2K  =  2.0,   DIS4K  =  0.2, 

ISPECL=  2,   DIS2L  =  2.0,   DIS4L  =  0.2, 

SMOO   =  1.00, 


VISCL  =  . F. , 


EPSE 

= 

0.35, 

$END 

$VISINP 

VISCJ 

= 

.F.,   VISCK  = 

NTURB 

= 

o, 

ITTYP 

= 

1/ 

ITDIR 

= 

3, 

JTLS 

= 

1, 

JTLE 

= 

42, 

KTLS 

= 

1, 

KTLE 

= 

45, 

LTLS 

= 

1, 

LTLE 

= 

50, 

TLPAR] 

0.3, 

$END 

$BCINP 
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NBC 

= 

8, 

IBTYP 

= 

1, 

14, 

12, 

54, 

12, 

32, 

30, 

40, 

IBDIR 

= 

3, 

2, 

-1, 

-2, 

1, 

-3, 

3, 

3, 

JBCS 

= 

1, 

1, 

42, 

I, 

1, 

1, 

38, 

38, 

JBCE 

= 

42, 

42, 

42, 

42, 

If 

42, 

42, 

42, 

KBCS 

= 

1, 

1, 

1, 

45, 

1, 

1, 

21, 

37, 

KBCE 

= 

45, 

1, 

45, 

45, 

45, 

45, 

26, 

40, 

LBCS 

= 

If 

1, 

1, 

1> 

1, 

50, 

1, 

1, 

LBCE 

= 

1, 

50, 

50, 

50, 

50, 

50, 

1, 

1, 

$END 

$SCEINP 

$END 
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APPENDIX  C.  FLOW  SOLUTIONS 
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Figure  C- 1 .  Mach  0.3  Inviscid  Convergence  History 
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Figure  C-2.  Mach  0.3  Viscous  Convergence  History 
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Figure  C-3.  Mach  0.4  Inviscid  Convergence  History 
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Figure  C-4.  Mach  0.4  Viscous  Convergence  History 
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Figure  C-5.  Mach  0.5  Inviscid  Convergence  History 
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Figure  C-6.  Mach  0.5  Viscous  Convergence  History 
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Figure  C-7.  Mach  0.6  Inviscid  Convergence  History 
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Figure  C-8.  Mach  0.6  Viscous  Convergence  History 
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Figure  C-9.  Mach  0.7  Inviscid  Convergence  History 
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Figure  C-10.  Mach  0.7  Viscous  Convergence  History 
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Figure  C-l  1.  Mach  0.8  Inviscid  Convergence  History 
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Figure  C-l 2.  Mach  0.8  Viscous  Convergence  History 
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Figure  C-13.  Mach  0.9  Inviscid  Convergence  History 
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Figure  C-14.  Mach  0.9  Viscous  Convergence  History 
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Figure  C-15.  Mach  1.2  Inviscid  Convergence  History 


xlO 


Figure  C-16.  Mach  1.2  Viscous  Convergence  History 
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Figure  C-17.  Mach  2.0  Inviscid  Convergence  History 
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Figure  C-18.  Mach  2.0  Viscous  Convergence  History 


91 


Figure  C-19.  Mach  4.0  Inviscid  Convergence  History 
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Figure  C-20.  Mach  4.0  Viscous  Convergence  History 
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Figure  C-21.  Mach  6.0  Inviscid  Convergence  History 
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Figure  C-22.  Mach  6.0  Viscous  Convergence  History 
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Figure  C-23.  Mach  6.0  Inviscid  Engine  Flow-Through  Convergence  History 
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Figure  C-24.  Force  Coefficients  vs.  Iterations  Mach  0.3  Inviscid 
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Figure  C-25.  Force  Coefficients  vs.  Iterations  Mach  0.3  Viscous 
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Figure  C-26.  Force  Coefficients  vs.  Iterations  Mach  0.4  Inviscid 
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Figure  C-27.  Force  Coefficients  vs.  Iterations  Mach  0.4  Viscous 
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Figure  C-28.  Force  Coefficients  vs.  Iterations  Mach  0.5  Inviscid 
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Figure  C-29.  Force  Coefficients  vs.  Iterations  Mach  0.5  Viscous 
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Figure  C-30.  Force  Coefficients  vs.  Iterations  Mach  0.6  Inviscid 
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Figure  C-31.  Force  Coefficients  vs.  Iterations  Mach  0.6  Viscous 
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Figure  C-32.  Force  Coefficients  vs.  Iterations  Mach  0.7  Inviscid 
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Figure  C-33.  Force  Coefficients  vs.  Iterations  Mach  0.7  Viscous 
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Figure  C-34.  Force  Coefficients  vs.  Iterations  Mach  0.8  Inviscid 
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Figure  C-35.  Force  Coefficients  vs.  Iterations  Mach  0.8  Viscous 
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Figure  C-36.  Force  Coefficients  vs.  Iterations  Mach  0.9  Inviscid 
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Figure  C-37.  Force  Coefficients  vs.  Iterations  Mach  0.9  Viscous 
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Figure  C-38.  Force  Coefficients  vs.  Iterations  Mach  1.2  Inviscid 
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Figure  C-39.  Force  Coefficients  vs.  Iterations  Mach  1 .2  Viscous 
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Figure  C-41.  Force  Coefficients  vs.  Iterations  mach  2.0  Viscous 
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Figure  C-40.  Force  Coefficients  vs.  Iterations  Mach  2.0  Inviscid 
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Figure  C-42.  Force  Coefficients  vs.  Iterations  Mach  4.0  Inviscid 
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Figure  C-43.  Force  Coefficients  vs.  Iterations  Mach  4.0  Viscous 
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Figure  C-44.  Force  Coefficients  vs.  Iterations  Mach  6.0  Inviscid 
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Figure  C-45.  Force  Coefficients  vs.  Iterations  Mach  6.0  Viscous 
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Figure  C-46.  Force  Coefficients  vs.  Iterations  Mach  6.0  Engine  Flow-Through 
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Figure  C-47.  Mach  0.3  Inviscid  Normalized  Surface  Pressure  and  Mach  Contours 


Figure  C-48.  Mach  0.3  Viscous  Normalized  Surface  Pressure  and  Mach  Contours 
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Figure  C-49.  Mach  0.4  Inviscid  Normalized  Surface  Pressure  and  Mach  Contours 


Figure  C-50.  Mach  0.4  Viscous  Normalized  Surface  Pressure  and  Mach  Contours 
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Figure  C-51.  Mach  0.5  Inviscid  Normalized  Surface  Pressure  and  Mach  Contours 


Figure  C-52.  Mach  0.5  Viscous  Normalized  Surface  Pressure  and  Mach  Contours 
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Figure  C-53.  Mach  0.6  Inviscid  Normalized  Surface  Pressure  and  Mach  Contours 


Figure  C-54.  Mach  0.6  Viscous  Normalized  Surface  Pressure  and  Mach  Contours 
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Figure  C-55.  Mach  0.7  Inviscid  Normalized  Surface  Pressure  and  Mach  Contours 
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Figure  C-56.  Mach  0.7  Viscous  Normalized  Surface  Pressure  and  Mach  Contours 
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Figure  C-57.  Mach  0.8  Inviscid  Normalized  Surface  Pressure  and  Mach  Contours 
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Figure  C-58.  Mach  0.8  Viscous  Normalized  Surface  Pressure  and  Mach  Contours 
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Figure  C-59.  Mach  0.9  Inviscid  Normalized  Surface  Pressure  and  Mach  Contours 


Figure  C-60.  Mach  0.9  Viscous  Normalized  Surface  Pressure  and  Mach  Contours 
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Figure  C-61.  Mach  1.2  Inviscid  Normalized  Pressure  and  Mach  Contours 
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Figure  C-62.  Mach  1.2  Viscous  Normalized  Surface  Pressure  and  Mach  Contours 
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Figure  C-63.  Mach  2.0  Inviscid  Normalized  Surface  Pressure  and  Mach  Contours 


Figure  C-64.  Mach  2.0  Viscous  Normalized  Surface  Pressure  and  Mach  Contours 
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Figure  C-65.  Mach  4.0  Inviscid  Normalized  Surface  Pressure  and  Mach  Contours 


Figure  C-66.  Mach  4.0  Viscous  Normalized  Surface  Pressure  and  Mach  Contours 
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Figure  C-67.  Mach  6.0  Inviscid  Normalized  Surface  Pressure  and  Mach  Contours 


Figure  C-68.  Mach  6.0  Viscous  Normalized  Surface  Pressure  and  Mach  Contours 
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Figure  C-69.  Mach  6.0  Engine  Flow-Through  Normalized  Surface  Pressure  and  Mach 

Contours 
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